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Abstract: 

Soft or collinear photon emission potentially poses numerical problems in the phase- 
space integration of radiative processes. In this paper, a general subtraction formalism is 
presented that removes such singularities from the integrand of the numerical integration 
and adds back the analytically integrated contributions that have been subtracted. The 
method is a generalization of the dipole formalism of Catani and Seymour, which was 
formulated for NLO QCD processes with massless unpolarized particles. The presented 
formalism allows for arbitrary mass and helicity configurations in processes with charged 
fermions and any other neutral particles. Particular attention is paid to the limit of small 
fermion masses, in which collinear singularities cause potentially large corrections. The 
actual application and the efficiency of the formalism are demonstrated by the discussion 
of photonic corrections to the processes 77 — > tt(7), e~7 — > e~7(7), and /i^/i^ — z/cZ/c(7). 
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1 Introduction 

Precision experiments with e^ beams, such as at LEP, at the SLC, or at future hnear 
coUiders, allow for an investigation of electroweak processes with a typical accuracy of 
some per cent down to some fractions of a per cent. An adequate description of such 
reactions — and a theoretical understanding of them that goes beyond a qualitative level — 
forces us to control higher-order corrections in perturbative predictions. An important 
source of such radiative corrections is due to the virtual exchange and the real emission of 
photons, or of gluons if quarks are involved. Although photonic corrections are formally of 
0{a) relative to the lowest order, leading to the naive expectation of ~1% as the typical 
size, the actual effects very often amount to ~10% or more. Apart from large kinematical 
effects caused by real photon radiation in particular processes, this enhancement mainly 
originates from coUinear photon emission off highly relativistic particles, such as e^ at the 
GeV scale, and from the corresponding virtual photon exchange. For initial-state radiation 
off electrons, this kind of correction is proportional to aln{'me/Q), where Q ^ me is a 
typical energy scale of the process. The remaining 0{a) corrections amount to one to 
a few per cent and have to be included in precision calculations as well. For per-cent 
accuracy even the leading 0{a^) corrections, or higher, can be relevant. 

In this paper we focus on the calculation of the full 0{a) correction that is induced 
by real photon radiation. Such calculations will be performed for practically all realistic 
obscrvables numerically, owing to the complexity of the squared amplitudes of the most 
interesting processes and the necessity of phase-space cuts. Usually the integration over 
the multidimensional phase space is performed by Monte Carlo integration. Thus, a 
linear increase in accuracy is roughly accompanied by a quadratic increase in the CPU 
time needed for the evaluation. In this context, the singularities of a squared amplitude 
cause problems. For example, the integrand of a bremsstrahlung process blows up if 
the photon energy becomes small, leading to the well-known logarithmic IR singularity 
in the phase-space integral. Following a frequently applied standard procedure, known 
as phase-space slicing, one introduces a small cutoff energy AE and integrates over the 
photon energy only down to AE numerically. The soft-photon part, E^ < AE, is known 
to factorize from the Born cross section, and the corresponding correction factor, which 
contains the IR singularity, can be calculated analytically. Since the results obtained 
this way are correct up to 0{AE/Q), precise predictions require rather small values of 
AE. For AE — > the numerical integration result grows like aln{AE/Q). Consequently, 
more and more CPU time is wasted in the precise calculation of this known singular term 
that cancels in the final result anyhow. Therefore, procedures that avoid such singular 
numerical integrations are desirable. 

Similar problems arise from collinear photon emission off a charged particle with mass 
m -C Q. Integrating over small emission angles 6 results in mass-singular corrections pro- 
portional to a\n{m/Q). Applying phase-space slicing, the coUinearity region is excluded 
by a small cutoff angle so that the singularity appears as a Q;ln(A^) contribution to 
the numerical integration result. The missing contribution from the region 6 < AO is 
related to the lowest-order cross section and can be obtained without singular numerical 
integration, similar to the IR case. Concerning the precision of the integration procedure. 
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plays a similar role as AE above, and a procedure that avoids the singular integration 
is preferable. 

Singular numerical integrations are absent in so-called subtraction methods. The idea 
of such methods is to subtract and to add a simple auxiliary function from the singular 
integrand. This auxiliary function has to be chosen in such a way that it cancels all 
singularities of the original integrand so that the phase-space integration of the difference 
can be performed numerically, even over the singular regions of the original integrand. 
In this difference the original matrix element can be evaluated without regulators for IR 
or collinear singularities, i.e. it is possible to apply powerful spinor techniques (see e.g. 
Ref. [ |l], 0, H and references therein) that have been developed for four space-time dimen- 
sions. The auxiliary function has to be simple enough so that it can be integrated over 
the singular regions analytically, when the subtracted contribution is added again. This 
part contains the singular contributions and requires regulators. In general, the statistical 
uncertainty of the finally obtained correction is smaller than the one of the correspond- 
ing result of phase-space slicing, because the absolute value of the numerical integral is 
usually much smaller for the subtraction method, owing to the absence of singular contri- 
butions. Unfortunately, the above requirements set highly non-trivial conditions on the 
subtraction functions, rendering the construction of a general subtraction procedure dif- 
ficult. Although various subtraction formalisms have been described for NLO corrections 
in massless QCD [ |^, to the best of our knowledge, up to now no general subtrac- 
tion method has been presented that is able to deal with massive particles in any given 
process. For the special case of heavy-quark correlations in hadron-hadron collisions, a 
subtraction procedure has been described in Ref. [ 

In the following we work out a rather general subtraction method for the treatment 
of photon radiation for any given process involving massive or massless, polarized or un- 
polarized fermions and any kind of neutral bosons. The inclusion of charged bosons is 
completely straightforward. Our method follows the guideline provided by the dipole for- 
malism, which has been presented by Catani and Seymour [ |^ for QCD with massless, 
unpolarized partons. Since the colour flow in QCD processes is more involved than the 
charge flow in electroweak processes, our presentation is simpler than the one in Ref. [ |^ 
in this respect. However, the generalization of the dipole formalism to arbitrary masses 
turns out to be highly non-trivial. Even in the limit of small fermion masses, which is of 
particular interest, there is an important difference between our approach and the sub- 
traction procedures of Refs. [ ^ Q for massless QCD partons. We consistently regularize 
IR and collinear singularities with finite masses, as it is commonly applied to photon 
radiation in electroweak processes, whereas the above-mentioned QCD studies are carried 
out in dimensional regularization. 

Although we treat only photon radiation explicitly, one should realize that the pre- 
sented results can also be used for gluon radiation in processes that involve only massive 
quarks as QCD partons; in this case the colour flow has to be handled as described in 
Ref. [ . Our work also represents a first step towards the generalization of the dipole 
formalism in QCD to include massive partons. 

The paper is organized as follows: in the next section we review the general structure 
of IR and collinear singularities, and describe the strategy of the subtraction procedure. 
In Section ^ we anticipate our results on the subtraction function and its integrated coun- 
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terpart for the special case of light fermions, in order to illustrate the structure of the 
formalism. The general results for arbitrary fermion masses are given in Section ^ where 
the details of the derivation are described, too. Section |^ contains the numerical exam- 
ples, including discussions of the photonic 0{a) corrections to the processes 77 — tt(7), 
e~7 — > e~7(7), and — > Uei^eil)- In Section | we discuss salient features of sub- 
traction formalisms and of the dipole approach. The discussion, in particular, includes 
comments on the implementation of phase-space cuts, some practical advice, and remarks 
on the partial generalization to QCD. Our summary is presented in Section ^ In the ap- 
pendix we provide important special cases, further details of the calculation, and the 
virtual photonic corrections to — > t'e^e- 

2 General strategy 

2.1 Preliminary remarks and conventions 

We consider photon emission in processes that involve arbitrary fermions and any 
massive neutral bosons. The initial state may also contain photons. The presented 
method remains applicable to reactions with more than one photon in the final state if 
only a single photon can become soft or collinear with a light fermion in phase space. Note 
that situations with more than one photon being soft or collinear correspond to corrections 
of or higher, relative to the lowest-order process without photon emission. In other 

words, the method to be described covers all kinds of real-photonic 0{a) corrections to 
processes involving charged fermions and any neutral particles. 

The relative charge and the mass of a fermion / are denoted by Q/ and ruf, the 
momentum and the helicity of / are assigned to pf and kj, respectively. Instead of the 
general indices /, /' for any fermions, we use the indices a, b only for initial-state fermions 
and i,j only for final-state fermions. Moreover, we define the sign factors ctj = ±1 for the 
charge flow related to the fermion /; specifically, we set ct/ = +1 for incoming fermions and 
outgoing anti- fermions, and af = —1 for outgoing fermions and incoming anti- fermions. 
Consequently, charge conservation of the whole reaction implies 

EQ/^/ = 0- (2.1) 

/ 

Since IR and collinear divergences are regularized by particle masses, we consistently 
work within four space-time dimensions. The invariant phase-space measure is abbrevi- 
ated by 

{2n)U^'^(^K-j2ki^. (2.2) 

In the following, A^i is the transition matrix element of the considered process that 
involves an outgoing photon with momentum k. The matrix element of the corresponding 
process without photon emission is denoted by A^o- For brevity, we explicitly write down 
only those momenta and helicities as arguments of Aii and Aio that are important in the 
considered equation. The collections of all momenta of the corresponding reactions are 
abbreviated by $1 and $0, and the respective phase-space measures by d$i and d$o- 



d(j){ki, . ..,kn;K) 



n 

.1=1 



(2vr) 



m, 
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2.2 IR and collinear singularities 

If the momentum k of the radiated photon becomes soft [k —>■ 0), the squared matrix 
element l-^iP? summed over all photon polarizations A^, becomes IR-singular and 
asymptotically proportional to |A^oP in the well-known form (see e.g. Ref. [ ^) 

EIA^.I\^o-i:«/-/W(J|^|A^oP. (2.3) 

where the sums on the r.h.s. run over all charged fermions of the reaction, and e is the 
positron charge. Eq. (|2.3|) is valid for all polarization configurations separately. The phase- 
space integral of ( |2.3| ) over the soft- photon region is logarithmically divergent. We choose 
an infinitesimal photon mass as regulator, yielding singular contributions proportional 
to a;ln(m-y) to the real 0{a) corrections. According to the Bloch-Nordsieck theorem [ 
^, these singular contributions cancel against IR-singular counterparts in the virtual 
corrections. 

Another type of singularity occurs in the limit of a vanishing fermion mass, m/ — 0, if 
the region of collinear photon emission from / is included in the phase-space integration. 
The squared amplitude |A1ip develops poles in (pfk) — > 0, leading to logarithmic singu- 
larities in the phase-space integral. The asymptotic form of J2x-, in the coUinearity 
regions is related to the squared amplitude |A1oP and well known [0, |^. Distinguishing 
between photon emission from outgoing and incoming fermions, we have 

E \M^ip,, QygiT\p^, k)\Mo{p^ + k; TK,)\', 

E ^^)\' pj^o ^'e'^r (Pa, k)\Mo{x,Pa, rKa)\', (2.4) 

where the signs r = ± account for a possible spin flip of the considered fermion. Whenever 
r appears more than once in products, we assume summation over r = ±. Note that we 
take helicity eigenstates as polarization basis throughout. The functions gf^^^^^^ are given 
by 



gl°+ '{Pi,k) 
gl°-\pi,k) 
ga!i{Pa,k) 



(1 - Zi f 



gL-\phk), 



2{p^ky 

1 

Xa 



2n 



where Pffiy) is the usual splitting function. 



1 +y2 

Pff{y) = 7^- (2.6) 
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If one is only interested in unpolarized fermions, the summation of g^^^^^^^^\M.Q{TKf)\'^ 

over r in (|2.4| ) reduces to (5^1°+*^™^ + gf^^^^^^)\M.o{K,f)\'^ . In the following, all formulas 
are written in a form that facilitates this replacement. The variables Zi and Xa are the 
fractions of the fermion energies that are kept by the fermions after photon emission, 

z--^ (27) 

While final-state radiation does not change any momentum other than pi in the hard 
scattering process, initial-state radiation scales the fermion momentum pa down to 
XaPa ^ Pa ~ thereby reducing the centre-of-mass (CM) energy of the hard scattering 
process. Integrating the squared amplitudes over a collinearity region yields contributions 
proportional to aln(mj) to the real 0{a) corrections^. According to the Kinoshita-Lee- 
Nauenberg theorem [ |10|, the mass-singular corrections aln(mj), which originate from 
final-state radiation, cancel against their counterparts in the virtual corrections in the to- 
tal cross section. Mass singularities from final-state fermions can only survive in specific 
distributions, such as distributions of invariant-masses that are built of fermion momenta 
only, i.e. without taking into account photon recombination. For fermions in the initial 
state the sum of real and virtual corrections remains mass-singular, and the Q;ln(ma) 
terms are a potential source of large corrections. 

The asymptotic relations ( |2.3|) and (|2.4| ) both relate the full squared matrix element 
|7Vli($i)P for the radiative process to |A^o('^'o)p5 which corresponds to the process with- 
out photon emission. Note that the phase spaces on which these functions are defined are 
different. In order to guarantee that both sides of ( |2.3| ) and ( |2.4| ) are defined on the phase 



space spanned by the momenta $i, one has to specify an appropriate mapping from $i 
to $0 that respects all mass-shell relations. The definition of such mappings is of central 
importance in the construction of a subtraction function |7WsubP that is parametrized by 
$1 and has the same asymptotic behaviour as J2x-, in the singular limits. 



2.3 The dipole subtraction formalism 

Our final aim is to perform the phase-space integral of J2x^ l-^iPi which involves 
IR and collinear singularities, without carrying out singular numerical integrations. The 
basic idea in a subtraction method is to subtract and to add the integral of an appropriate 
subtraction function |A^subP, 



J d<l>i J2 l-^il' = / 1 [E l-^il' - l-^subl'] + / d$i |^,,b 



2 

5 



(2i 



where |A^subP possesses the same asymptotic behaviour as J2x^ l-^iP in the singular 
limits. Specifically, we demand 



^There are also corrections of the forms aln(m/) In(m^) and aln^(m/), which originate from soft 
photons. These corrections, however, always cancel against virtual corrections. 
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where i and a label all outgoing and incoming light fermions. Owing to (2^), the phase- 
space integration of the difference (J2\-y l-^iP" l-^subP) in (|2.8| ) is non-singular, i.e. it can 
be performed numerically without regulators. The singular contributions of the original 
integral / d$i J2x^ are completely contained in / d$i|A^subP- If |A^subP is chosen 

appropriately the singular integrations can be carried out analytically. To this end, the 
phase-space integral / d$i is factorized into a part / d$o connected to the non-radiative 
process and a part J[dk] connected to the photon phase space, 

d$i = j d%® j [dk]. (2.10) 

The sign "®" indicates that this factorization is not an ordinary product, but may contain 
also summations and convolutions. Since X^a^ |-A^i('^'i)P is related to |A^o('^'o)P of the 
non-radiative process in the singular limits, the subtraction function |A^sub('^'i)P can be 
defined in such a way that it depends on the momenta of $o only via |A^o(^'o)P- The 
integration variables of /[dfc] occur only in the remaining terms of the subtraction function. 
Since those terms are process- independent, the singular integration of |7Wsub('^'i)P over 
[dfc] can be performed analytically once and for all. Finally, the integral / d$i X^a^ 
takes the schematic form 

j d$i i-^ii' = / d-^"! [e i-^iT - i-^subp) + /d<io ® (1 m lA^subr) . (2.11) 

The integrations over d$i and d$o on the r.h.s. are free of singularities, and thus are 
well-suited for numerical evaluations. Since the singularities in /[dfc]|A^subP are controlled 
analytically, they can be easily combined with their counterparts in the virtual corrections. 

We have seen that the subtraction function |A^subP has to obey two non-trivial condi- 
tions. It must possess the asymptotic behaviour given in ( |2.9D , and it must still be simple 
enough so that it can be integrated over the singular regions analytically. Note that all 
the collinearity regions of phase space overlap and have the IR part {k — > 0) in common. 
Therefore, the naive sum of all coUinear singularities, which are proportional to Qj, leads 
to an overcounting of the IR singularity, and thus cannot be used in the subtraction func- 
tion. In the following we show how this overcounting is avoided and how the subtraction 
function is constructed within the dipole formalism. In contrast to Ref . [ |^ , where this 
formalism is described for massless, unpolarized partons in QCD, we have to take care of 
fermion masses and polarizations. 

The subtraction function |A^subP is constructed from auxiliary functions gfj/^l, which 
are labelled by pairs of different fermions / 7^ /': 

|-Msub($i)|' = - E Qf^fQi'^f^'ap^liPf^Pf^f') 1-^0 {^o,fr;rKf)\\ (2.12) 

Since only the kinematics of fermion / gives rise to singular contributions in the sub- 
traction function, / is called emitter, whereas /' is called spectator. The summation over 
r accounts for the helicity flip of the emitter /. The auxiliary functions gfj^l have to 
possess an appropriate asymptotic behaviour. In the IR limit one globally demands 



fc->0 



Pfk 



m 



Pfk+pfik Pfk, 
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(sub) _ (sub) _ 

'ij • 9ia 




Figure 1: Effective diagrams for the different emitter/spectator cases. 

gf;};l{pf,Pf,,k) ~0{1). (2.13) 

In the colhnear hmits one demands separate conditions for final- and initial-state fermions 
f = i,a: 

9i;^liPa,Pf', k) g^^ip,, k), (2.14) 

where /' can be outgoing or incoming. The analytical form of gfj^l is, of course, not 
uniquely determined by the asymptotic conditions. A convenient choice for these auxiliary 
functions, which are graphically represented by the effective diagrams of Fig. |l|, is given 
in the next sections. 

Finally, we have to specify the conditions on the momenta to be inserted in |A^oP in 
( p.l2| ). As explained above, it is necessary to define a mapping $o,//' from the momenta 
of $1 to the ones of $o that respects all mass-shell relations. The symbol $o,//' indicates 
that different mappings are used for different pairs //'. Denoting the momenta of / and 
/' in $0,//' by and pf, and writing kn and /c„ for the remaining momenta in $i and 
$0,//'; respectively, we require 

P/^P/, Pf'^oPf'^ ~^n^,kn (2.15) 

for the IR limit and 

p^p;^oP^ + ^' p-pj^o^-p-^ Pf'p;^oPf'^ ^«p;^o^- (2-16) 

for the colhnear limits. The mass-shell conditions p'j = p'j = rnj-, pj-, = p'j, = m^,, and 
kn = k'^ = have to be fulfilled for arbitrary photon momentum k. 
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Using ( p. 13 ) and ( p. 14 ), it is rather easy to check that the subtraction function ( 2.12|) 
possesses the asymptotic behaviour required in ( p.9|) . For the IR hmit, the asymptotic 
relation is verified upon inserting ( p.l3| ) into ( p.l2| ) and rearranging the terms in J^fj^f'- 



J2 <5/^/Q/'^/'— r 



m. 



Pfk + pfik pfk 



Pf'Pf 



.2 rn} 



fj' 



{Pf'k){pfk) 

Pf'Pf 
{Pf'k){Pfk)' 



f 



:2.i7) 



Note that charge conservation ( |2.1| ) was used in the form 

Yl Qf'^f = -Qf^f 
f'if'^f) 



:2.i8i 



in the term proportional to m?j. The arguments of |A^oP iii |-^subP behave in the desired 
way owing to ( p. 151 ). The asymptotic relation ( |2.9D for the coUinear limits follows after 
inserting the conditions ( |2.14]) into ( |2.12|) and using again charge conservation ( |2.18| ). The 
correct behaviour of the momenta of |A^oP in |-^subP is guaranteed by (p.l6|) . 



3 Subtraction functions and integrated counterparts — special case of light 
fermions 

Before we turn to the treatment of the general case of massive fermions in the next 
section, we first describe the dipole subtraction formalism for light fermions, i.e. we neglect 
fermion masses in this section whenever possible. In this way, the structure of the 
formalism becomes clear without being obscured by all kind of complications that are 
related to particle masses. Moreover, this section provides a condensed instruction to the 
formalism for light fermions, since details of the method that are only relevant for its 
derivation are also postponed to the next section. 

3.1 Final-state emitter and final-state spectator 

We define the auxiliary functions gij^r \ which correspond to a final-state emitter i 
and a final-state spectator j, by 



gl:]'^^\pi,Pj,k) = 0, 



where the variables i/ij and Zij are given by 

Pik 



Vij 



PiPj +Pik + pjk 



Zij {I 



y^j) 



- 1 - Zi 



z 



pipj 



PiPj +Pjk 



(3.1) 



(3.2) 



^The dipole subtraction formalism for this important special case has been worked out independently 
by M. Roth [ 0. Comparing both approaches, we find full consistency. 
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Since we assume rriij —>■ 0, the explicit mass terms in gl"^^^ and in glj^r^ are negligible 
in the difference (|A^subP — J2x^ l-^iP)- particular, this implies that g^f^^^ vanishes. 
Note, however, that those mass terms are relevant in the integration of glj^r^ over the 
photonic part of phase space (see next section). It is straightforward to check that the 
functions glj^^^ of obey the asymptotic conditions (|2.13|) and ( p.l4| ) in the IR and 
coUinear limits, in which we get 



,2 

with Zj from (12. 71). For the evaluation of 



-^0 {^o,ij) "we have to define the mapping l>o,ij 
from $1 to $0- Of course, it is desirable to leave as many momenta unchanged as possible. 
Therefore, we redefine only the momenta of / and /', and leave all other momenta 
unaffected, kn = kn- The momenta pi and pj are chosen as 

Vij H ~u 1 

—Pj, Pj = -— 

Vij ^ Vij 



P^ = P^ + k^- T^Pi ' Pi = 1 P'r (3-4) 



The on-shell relations = = and the validity of the required asymptotic behaviour 
( p.l5| ) and ( |2.16| ) can be checked easily. Moreover, momentum conservation. 



Pij = Pi+Pj + k = Pi+ Pj, (3.5) 

is trivially fulfilled, and P^j > holds for all phase-space points. The above defini- 
tions comprise all ingredients for the evaluation of the ij contribution to the difference 
(J2x^ 1-^1 P — l-'^subP)? which is integrated over the full phase space numerically. We 
recall that this integration can be performed with vanishing photon and fermion masses. 

The construction of the above contribution to the subtraction function entirely fol- 
lows the pattern of Ref. [ H for massless QCD partons. The same applies to the other 
emitter/spectator cases. Note that we have added a factor of 1/(1 — yij) in g^j^^^ that is 
not included in the approach of Ref. [ ^. This factor, which is introduced for convenience 
(see massive case in Section [4.1|) , affects only non-singular contributions. 

The differences between Ref. [ ^ and our approach for light fermions become apparent 
in the analytical integration of the subtraction function over the photonic part of phase 
space, where both IR and coUinear regulators have to be taken into account. The photonic 
part of phase space is defined by extracting the phase-space measure d$o,ii! which is 
spanned by the momenta pi, pj, and from the full phase-space measure d$i, which 
is spanned by the momenta Pi, pj, k, and The details of this splitting and of the 
integration over the photonic part can be found in Section ET] and in App. O. Denoting 



the integral of glj^r^ over the photonic phase space by G'i^^^\ and including an appropriate 
normalization factor, the ij contribution \ Aisuh,ij{^i)\^ to the phase-space integral of the 
subtraction function reads 

J d$i |A<sub,i,($i)r = -^Q^a^QJaJ J d%,,Gi;^^\P^^)\Mo{pi,pj;rK,)\', (3.6) 

where a = e^/ (Ait) is the fine-structure constant. The functions G^^^^^ are explicitly given 
by 

= ^(4'^?) - y + 1' G^{P^ = (3.7) 
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where the singular terms are contained in the function 



.,P^™^,^>.<9!)>.<9i)..(9i)-i...(;!).i.(9!)^ 



(3i 



As required, only the emitter mass gives rise to logarithmic singularities, whereas the 
spectator mass rrij can be set to zero exactly. The spin-flip part G[^^!^^ is non-vanishing 
and entirely induced by photons that are emitted coUinearly to the emitter i. Since all 
singular terms are factorized into gI^^^\ the masses m-y, rrii, and rrij can be set to zero 

everywhere in ( |3.6| ) apart from G\^^^\ 

3.2 Final-state emitter and initial-state spectator, and vice versa 

Emitter/spectator pairs from the final/initial state and vice versa, i.e. the cases ia and 
ai, always occur in pairs for a given process. Since the kinematics is identical in both 
cases, we treat them in one go. The corresponding auxiliary functions gff^l are given by 



yia,+ 



„{sub) 



(sub) 



1 

{Pak)x 



2 



2 ^ia ^it 



{Pi,Pa,k) = gaC^\pa,Puk) = 0, 



(3.9) 



with the variables 



PaPi+Pak-pik 



PaPi 



(3.10) 



PaPi+Pak PaPi+Pak 

The desired asymptotic behaviour (p.l3|) and ( p.l4| ) in the singular limits, which imply 



(3.11) 



can be verified easily. The modified momenta Pi and Pa of the sets $o,ja = ^o,ai are chosen 
as 



Pf=pf + k>'-{l- X,a)p^„ = X^aP^a, 



(3.12) 



and the remaining momenta kn coincide with the corresponding momenta kn of $i. The 
on-shell relations pf = Pi = 0, the required asymptotic behaviour ( [^.15| ) and ( ^.16| ), as 
well as momentum conservation, 



Pi + k - Pa= Pi- Pa 



(3.13) 



can be checked easily. For massless fermions i and a we always have P^^ < 0. This 
completes the definitions of all quantities for the evaluation of the ia and ai parts of 
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the difference — l^subP); which can be performed numerically with vanishing 

photon and fermion masses. 

The analytical integration of the ia and ai parts of the subtraction function is more 
involved than in the ij case, since the modified momenta pi and pa correspond to a new 
initial state. In the following, we consider a scattering reaction with the two incoming 
light-like momenta Pa and pf,. Owing to Pa — XiaPa, the CM frames of Pa +Pb and Pa +Pb 
arc related by a boost along the beam axis. The strength of this boost is determined by 
Xia, which is the ratio of the corresponding squared CM energies s and s, 



S = 2{paPb) = 'iXiaipaPb) = XiaS. 



(3.14) 



The photonic part of the phase space, which results from the extraction of the phase- 
space measure d$o,m from the full measure d$i, involves an integration over Xia- This 
integration over Xia cannot be carried out analytically, since the complete phase space 
spanned by the new momenta $o,ta implicitly depends on Xia via the CM energy a/I. Thus, 
the integral over the photonic part of phase space is written in terms of a convolution 
over Xia- Including an appropriate normalization, the contributions to the phase-space 
integral of the subtraction function read 



a 



Zn 



X 



/ dx d^o,iai^) - Qff^liPia, ^) Mo(xpa + Pia,XPa]TKf 

«y J t£ 



(3.15) 



with //' = ia, ai for the two different cases. In this convolution x plays the role of 
Xiai and the argument of the phase-space measure d$o,ia(3^) indicates that each value of x 
determines a different phase space. The momenta to be inserted in | AloP are pi — xpa+Pia 
and pa = xpa, where Pa is fixed, and Pia varies with the phase-space point in ^o^ia{x). 

Since the distributions Gff^]-{Pij^,x) become IR-singular at the point x ^ 1, the 
convolution is not yet suited for a numerical evaluation. A possible way out is provided 
by the application of the [...]+ prescription to this distribution. 



''dx[f{x)Ug{x) = ['dxf{xMx)-g{l)]. 

') Jo 

Using this trick, the IR singularities in the endpoint contributions 

Gf?:i{Pi)^£dxgp}:i{PLx) 

are separated from Q^f^!^l{P^g^, x) , and the convolution reads 



(3.16) 



(3.17) 



d$l |Alsub,//'(^l)P = -^Qa<^aQi(Ti 




d^O,iaiPla,x) 



X 



^(sub) 



(PL^) 



ld^oAPlA)Gff^liPL^) 



Mo(xpa + Piai XPa] TKf) 
Mo(pa + Pia,Pa]rKf^ 



+ f d^o,a{PL 1) G(7,'J(i^2j \Mo{pa + P^a,Pa, TKf) 



(3.18) 
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Note that we have included P^^ as additional argument in the phase-space measure 
d^Q ia{Pia, x) , in order to signalize that we have kept fixed during the integration 

7(sub) x-y(sub) 



TBI) all singular 



over X in the calculation of the endpoint contributions G^^j"l{Pl^] 

contributions are factorized into Off/^l or G^jj}'l so that this equation is well-suited for 



numerical evaluations, and the photon and fermion masses can be set to zero everywhere 



apart from Qff^l and G 



^(sub) 



Finally, we give the explicit form of the distributions Q^j^, l iP^ 



7 (sub) 



7 (sub) 



(P2 

V It 



(P^ 

\ IC 



,X 



X 



X 



2 In 



In 



X 



1 — X 
Pi 



mix 



1 — X 



ln(2 



x] 



'(sub). 



P^ x) 



J (sub) 



(P2 

V u 



X 



0, 

1 



and the corresponding endpoint parts G^^j}'^^{Plj) , 



71 

T 



[l + x) ln(l 



X) 



(3.19) 



G^{P^ = C(]P^,n^ + --^ 



TT 



^(sub)/p:^N _ ^^SUDj/p 



(sub) , 



ia) 



1 

2' 



(3.20) 



where C is the function defined in ( p.8|) , which contains the logarithmic singularities. Only 
the emitter masses lead to mass singularities, as it should be. From the general discussion 
of mass singularities in Section |2.2| , it is also clear that the singularity of final-state emitter 
appears only in the endpoint contribution Gff^. For an initial-state emitter also the 

distribution Q^ai^^ contains a mass-singular part, which is proportional to the splitting 
function Pff{x). The spin-fiip contributions are regular, and the one for a final-state 
emitter is completely contained in the endpoint part. 



3.3 Initial-state emitter and initial-state spectator 

For an emitter a and a spectator b from the initial state we introduce the variables 

PaPb-Pak-pbk Pak 



Xab 



PaPb 

and define the auxiliary functions ga^r^ by 



Vab 



PaPb 



(3.21) 



9ab!+{Pa,Pb,k) 



{Pak)Xab 



1 - Xab 



1 — X, 



ab 



9ab^^\Pa,Pb,k) = 0. 



(3.22) 
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They posses the required asymptotic behaviour in the singular hmits, which are charac- 
terized by 

^ab ^ 1, Vab ^ 0, Xab ^J^q Xa, Vab ^J^q 0. (3.23) 

Following the guideline of Ref. [ the construction of the modified momenta $o,afe, which 

are used to evaluate A4o ($o,ab)| in ( p. 121) , differs from the previous cases. Instead of 
changing only the emitter and spectator momenta, we now keep the spectator momentum 
Pb fixed and change all outgoing momenta kj other than k. Note that kj also includes the 
momenta of neutral outgoing particles, i.e. we have 

Pab = Pa + Pb-k = J2 kj ■ (3.24) 

j 

The new momenta 

P^,=XabP'a, P!:b = XabP'a+P'b (3-25) 

are chosen in such a way that p'^ = and P^j, = P^^. In the IR limit they obviously tend 
to Pa and Pab, respectively; in the coUinear limit approaches XaPa with Xa from (p. 7]). 
The individual momenta kj are modified by a Lorentz transformation as follows, 

= Ki^Xj (3.26) 

with 

A/i _ /I (-Pgfe + Pab)^iPab + Pab)u '^PgbPab^u /„ r,«\ 

^^u-9u p2 , p p + p2 ' y'^-^') 

J^ab + J^abJ^ab ^ab 

SO that the mass-shell relations = are retained. The necessary condition K^^K^f = 
Qj^p and the relation kj = Pat are easily checked by direct calculation using P^^ = P^^,. 
The above relations comprise the necessary input for the construction of the ab part 
|-^sub,abP of the subtraction function. 

Concerning the analytical integration of the ab contribution to the subtraction function 
over the photonic part of phase space, the situation is similar to the previous section. The 
incoming momenta Pa and pb are modified in an analogous way, namely Pa is scaled down 
to pa by the variable Xab, while pb is kept fixed. The CM frames of pa +Pb and pa +Pb are 
related by a boost along the beam axis, and the squared CM energies s and s are related 
by 

S = 2{paPb) = 2Xab{PaPb) = ^abS. (3.28) 

Note also that P^j, = s = XabS, owing to definition ( p.21|) . The separation of the photonic 
part of phase space again leads to a convolution over x = Xab for the integrated ab 
contribution |A^sub,afeP to the subtraction function. 



Jd^l |A^sub,a6($l)P = -^Qa(^aQb(^b 

X £ dx^gl;','^i\s,x) I d^oA^) \Mo{xpa,kn{x);TKa)\\ (3.29) 



The argument of the modified momenta fc„(x) indicates that the new phase space implic- 
itly depends on x. For the numerical evaluation of the convolution, it is appropriate to 
separate the IR-singular endpoint part 



f'dxg;^,:i\s,x) (3.30) 
Jo 
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from the distribution Gai^r^ with the help of the [...]+ prescription. The numerically 
accessible form of the convolution is 



d$l |A^sub,afe($l 



a 
"2^ 



QaO'aQbC^b 



where the distributions are given by 



1 /" ~ ~ 
- J d$o,afe(s,a;) Mo[xpa, kn{x);TKa 

d<lo,ab(s, 1) Mo{pa,kn{l);TKa 
) / d<lo,afe(s, 1) Mo{pa,knil);rKa) 



GSI\s,x) 



P 



ff 



x] 



In 



(4)- 



and the endpoint parts read 



(sub) , 
ab,+ \ 



0^:-\s,x) 



1 — x, 



(3.31) 



(3.32) 



(3.33) 



Note that the original squared CM energy s is kept fixed in the integration over x that 
defines the endpoint contributions in ( p.30| ). This is also indicated in the phase-space 
measure d^o,ab{s, x), which is to be parametrized for fixed s and x. The mass singularities 
of the emitter a are completely factorized into 0^1^^ and G^^^^ so that the convolution 
over X can be carried out numerically for vanishing photon and fermion masses. Of course, 
the spectator mass is set to zero everywhere. 



4 Subtraction functions and integrated counterparts — general case 

In this section we turn to the case of arbitrary finite fermion masses. Here we include 
also details of the derivation, which have been omitted in the previous section for brevity. 
The anticipated results for light fermions can be obtained from the general ones of this 
section by carefully expanding the corresponding formulas for small fermion masses. 

Moreover, it is phenomenologically important to consider the case of light fermions 
only in the initial state, which is of particular interest for e"'"e~ collisions at high energies, 
as observed at LEP or the SLC In this case, the dipole formalism is also considerably 
simpler than for general fermion masses. The results of the corresponding expansion are 
listed in App. [A.l . 

4.1 Final-state emitter and final-state spectator 

In order to construct the contribution to |A^subP corresponding to an emitter i and a 
spectator j from the final state, we have to define the auxiliary functions gf^^\pi,pj, k) 
and the embedding of the momenta pi and pj into the phase space spanned by pj, pj, and k. 
This embedding has to respect the mass-shell conditions pf = pf = mf, p'^ = p'^ = m'j, and 
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k'^ = m?, where the photon mass is taken to be infinitesimal in the final result. In the 



following we make use of the variables yij and Zij of (|3.2| ) and introduce the abbreviations 
P^,=Pi+p, + k, P^^ = P^^ ~ ml - m] - ni^, A,, = A(4, m^, m^), (4.1) 
where 



define the auxiliary functions 

Rij{y) 



A(x, y, z) = + + z- 2xy - 2xz - 2yz. (4.2) 

^2 p2 \ 
ij-i ij-i ■^ij 



In the physical phase space we always have P^^, Pf^, Xij > 0. For later convenience, we 



(2m| + /^2_^2^)2-4/^>2 
Ajj 

2m2(2mf + /^2) ^ 



r,(.)^ 1-^-3^^. (4.3) 

Their actual form is convention except for their behaviour near y = 0, where they are 
regular with Rij{0) = rij{0) = 1 for = 0. 

Using the above quantities, we define the functions g^T^^ as follows: 



^^'^^^P-P^^''^ (p.A;)i?,,(y,,)[l-^,,(l-y,,) ^ ''^ 
^ ,N rn] (1 - z^j^ r,j{y 



1 



Pik 



- 9iT~{Pi,Pj,k), 



It is straightforward to check that these functions obey the asymptotic conditions (|2.13|) 
and ( |2.14| ) in the IR and coUinear limits, in which yij and Zij behave as given in ( ^.3| ). Note 
that the limits (|3.3|) implicitly assume 0; the collinear limit additionally requires 

rrii 0. In the mapping $o,ij from $i to $o we leave all momenta kn other than pi, pj, 
and k unaffected, as in the case of light fermions. The momenta pi and pj are chosen as 



I } rV^i p2 Vj^ 9p2 ii' 

^A((p. + fc)2,/^2^m|)^ ^ 

pr = i^';-pT (4.5) 

The on-shell relations = and = mf directly follow by expanding the squared 
momenta, and momentum conservation Pij = pi+pj + k = Pi + Pj is fulfilled by definition. 
Moreover, the validity of the required asymptotic behaviour ( p.l5| ) and ( |2.16| ) is obvious. 

The relations given above include all ingredients that are necessary to calculate the 
ij contribution to the subtraction function |A1sub('^'i)P- The phase-space integral of the 
difference {J2x^ lAliP ~ l-'^subP) is non-singular in the IR limit, and also in the collinear 
limit, which occurs for rrii —>■ 0. Therefore, this integral can be evaluated with = 
everywhere, and with rrif = for light fermions. However, we need the dependence on 
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and light fermion masses mj in the integral of |A^subP over the photon phase space, 
which is calculated next. 

To this end, we separate the phase space spanned by the momenta pi and pj from the 
one that is spanned by pi, pj, and k as follows: 



d(j){pi,pj, k; Pij) = J d(f){pi,pj; Pij) J [dkiP^j, yij, Zi^ 



(4.6) 



Our aim is to perform the integration of the ij part of the subtraction function over the 
remaining phase-space variables contained in the measure [dk{P^j,yij, Zij)]. The explicit 
form of this measure, which is derived in App. |B|, is given by 



P.4 



4(27r)3,^/v-yo 



2n 



dipi 



y2 



yi 



dy^j (1 - yij) 



22{j/ij) 



dz. 



(4.7) 



The angle (pij is the azimuthal angle of Pj with respect to the Pj axis in the CM frame of 
Pij. The integration boundary for the variables yij and zij is given by 



2/1 



2/2 



P.2 



Zi,2{yi: 



{2m1 + Pfjyij){l - yij) T ^vlj - Vi^Kj Rjjjy 
2{l - y,.j){m} + m^^ + Pf^y,,) 



(4.8) 



Since the integrand gij'^r\piiPji ^) -^o {^oaj', Ti^i) of the phase-space integral does not 
depend on the angle yjjj, the integral over ^pij simply yields a trivial factor of 27r. Moreover, 

M.0 {^o,ij) is independent of yij and Zij so that the integrations over yij and % only 



concern the auxiliary functions giTr \ ^ind we define 



P.4 



y2 



(^Vij (1 - Vij) 



yi 



Z2(yij) 
21 fe) 



(4.9) 



While the integration over Zij is very simple, the one over yij is non-trivial, but can be 
performed analytically. Details of the calculation can be found in App. 0. We obtain 



^(sub) / p-2 



In 



mf 



-21n(l-a^) + ^ + - + 



PI 



In(ai) In 



+2Li2(ai) +4Li2 



TT 



— - AU2{-^/a^) + - In (ai) - — 



(sub) 



P2.) 



P.4 



In 



+ 



+ (1 



mj{Pij + 2m: 



P.4 



+ 



2a3 

[4arctan(a3) — vr] 

2m?mi 



+ 
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+ 



1 2 
miVTLj 



In 



'2mimja^^ 



(4.10) 



with the shorthands 



Pij + - ^Xij - ^Xij rrii 
^1 = 7^' ^2 = — 7=, as = -j=F , (4.11) 

and the usual dilogarithm Li2(x) = — Q dtln(l — xt)/t. The contribution of the ij part 
|-^sub,ii(*&i)P to the phase-space integral of the subtraction function is formally the same 
as given in ( |3.6| ) for light fermions, where we have to insert the functions Cf^^'' of ( |4.10| ) 
for finite fermion masses. Expanding these functions for rrij j —>■ 0, we obtain the results 
of 

4.2 Final-state emitter and initial-state spectator 

For a final-state emitter i and an initial-state spectator a we keep the definitions ( |3.10|) 
of the variables Xia, and introduce the abbreviations 

Pia=Pi + k-pa. Pl = Pl-ml-m1-ml^, Xia = X{PL ^1 ^^i) ■ (4.12) 

In the following, we only consider fixed momenta Pia that obey 

Pl<{ma-mif, (4.13) 

because other values of P^^ do not admit the limits k ^ or pik ^ 0, as can be checked 
easily. In particular, ( |4.13| ) ensures that P^^ < 0. Moreover, we introduce an auxiliary 
parameter xq with < < 1 that specifies the lower limit on Xia for which the subtraction 
function will be applied. We are forced to deviate from the simple choice Xq = 0, because 
Xia — is not allowed for all configurations of P^^, ma, and nii. More precisely, one has 
to require 

-P^ I — 

Xo> X = ^ if < JPfa <ma- TTli, (4.14) 



2ma [rria - \l Pta 

which is only possible for Pf^ > and rria > rrii. Otherwise we can take any value for xq 
with < Xq < 1. For instance, it is possible to set Xq = for vanishing fermion masses, as 
done in Section |3^ . For P^^ > {rria — ^iY ^ia < the subtraction functions g^^^^ are 
set to zero consistently. Note that both IR and collinear singularities appear at Xia —>■ 1, 
i.e. applying the subtraction function for xq < Xia < 1 correctly cancels these singularities. 
The final results on observables must not depend on xq, which will not be further specified 
in the following. Checking the xq independence of observables is a non-trivial check on 
the complete subtraction prodecure. 

It is convenient to introduce the auxiliary functions 



P . , J{Pl + 2mlxy~AmlPlx' 

V Xi 



, ^ 1 , ^a(^i + 2m2)l-a; ,^ ^ , 

Tiaix) = 1 + ^ , (4.15) 
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with -Ria(l) = fia{^) = 1 for vanishing photon mass m^. The subtraction functions for 
finite fermion masses read 



9iT-{Pi,Pa,k) 



{Pik)x 



2 



2{p,kf 



- 1 - Zi, 



{Xia) 



Pik 



-9, 



(sub) 



{Pi,Pa,k), 



(4.16) 



They possess the desired asymptotic behaviour ( p.l3| ) and ( p.l4| ) in the singular hmits. 
The behaviour of Xia and Zia for /c — » or pjA; is given in ( 3.11 ). The auxihary 

momenta $njV), which are needed to evaluate 



M ('^o,ja) I , are constructed in a similar 



way as for light fermions. The emitter and spectator momenta are given by 



P'a 



x{{Pi + ky,p, 

Ft lai 



L ^1 



P'a 



PiaPa 



PL\ + 



P2 



2R 



iai 



(4.17) 



and the remaining momenta k^ coincide with the corresponding momenta k^ of $i. The 



on-shell relations pf 



Pi 



mj and 



pI 



can be verified easily. Momentum 



conservation Pia = Pi + k—pa =Pi—Pa and the validity of the required asymptotics (|2.15|) 
and ( |2.16| ) are obvious. 

Using the above relations, the ia part of |A^sub('^'i)P can be evaluated. Concerning 



the role of the masses and m/ in the difference {J2xj l-^iP ~ |A^subP), the remarks 
of the previous section apply as well. 

As already explained in Section p.2| , the photonic part of phase space, which is ob- 
tained by separating d$o.ia from d$i, involves an integration over the variable Xia- This 
integration over Xia turns into a convolution over x = Xia in the integration of the subtrac- 
tion function, where x determines the CM energy of the reduced phase space spanned by 
the momenta $o,ia- For the explicit splitting of phase space, we decompose the momentum 
Pia = Pb — Kia iuto the secoud incoming momentum and the total momentum Kia of 
all outgoing particles other than i and the inspected photon. The phase-space separation 
is defined by 



j d(j){pi,k,Kia;Pa+Pb)0{xi. 



Xo) 



dx 



xo 



d(j)[pi{x),Kia;Pa{x) +Pb) / [dk{P^^,X,Zia)] 



(4.18) 



with Xi given below. The x-dependent momenta 

1 



Pai^) 



pf{x) 



Pia (•^) 



Pai^) 



„ Pi + 2mlx 
xpf^ + _ ^„ " PH 



2Pi 



Pfa + K- ^1 pM 

2Pl 



P^ 



(4.19) 



result from pa and pi upon substituting (paPia) ~^a ~ PiaK'^^)^ which eliminates 
iPaPia) in favour of x and Pf^. Note that the substitution also concerns {pi + k)"^ = 
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{pa + Pia)"^ in ( |4.17| ). If we now try to reconstruct p^ix) from x and Pa, as it was possible 
for nia = in Section we find that the knowledge of x and Pa is not sufficient for 
TUa 7^ 0. This complication is due to the fact that the boost relating the CM frames of 
Pa+Pb and Pa+Pb does not simply go along the beam axis for finite rria- We will come back 
to this problem and to its solution at the end of this section and proceed by performing 
the integral of the contribution to the subtraction function over [dfc(Pj^, x, Zia)]- In App. ^ 
the explicit form of this measure is derived. The result is 



Pi, 



4(27r) 



X^ 



z\{x) 



2tt 



dip 
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(4.20) 



where (p^ is the azimuthal angle of the photon in the CM frame oi pi + k. The function 

Pia{s) reads 



\ A(s,m2,m2)' 



(4.21) 



where s and s denote the squared CM energies of Pa + Pb and Pa + Pb, respectively. The 
integration boundary for Zia is given by 



Zi^2{x) 



^P^a - x{Pl - ml)] 
From this relation we read off that the maximal value of x is given by 



Xl 



P2 



2mjm^ 



1 - + 0(m2 



\P^ 



1' 



(4.22) 



(4.23) 



While the integration of gl^^\pi, Pa, k) Aio ("^0,^; T/ti) over ip^ yields a trivial factor of 



2tt, the integration over Zia depends on the actual form of g 



(sub) 



Defining 



P4 



dzia gff}^l{pf,pf',k) 



for //' = ia, we obtain 
Pi 



7 (sub) 



(PLx) 



2^ \iaRia (•^m) '^^l(^ia) 

2 — X — Zi{x) 



(4.24) 



7 (sub) 



{PLx) 



m, 



\a Ria{x){l — x) 
+ [Zi{x) - Z2{x) 

xria{x) 



2 In 



1 + 



2 — X — Z2{x) 
Zi{x) + Z2ix) 



2m?x 



Plil-x) 



ot?iPLx 



/A~(l - x^Riaix) 

+ [Zi{x) - Z2{x)] 



jln 


Z2{x) 




_zi{x)_ 



Zi{x) + Z2{x) 



(4.25) 



after a simple integration. The function gIo^^^ is singular in the limit x —>■ 1 for = 0. 
Since our aim is to perform the convolution over x numerically with = 0, we separate 



19 



the singularity at x — 1 by introducing the [...]+ prescription. Considering that our 
lower integration limit Xq can be different from zero, we write 



(sub) , 



+ 



^(sub) / p2 



X 



with 



(4.26) 



(4.27) 



ia. In writing ( ^.26| ), we have already used that the photon mass 

(sub) 



where //' 

can be set to zero, and thus Xi set to 1, everywhere apart from G'f^^\ which contains 
the IR singularity. The integration over the distribution Q^^^\p^^,x)6{x — Xq) ^ can be 
performed with = and, if desired, with m/ = for the fermion masses. The endpoint 
contributions Gi^^^\Pl^, Xq) are obtained after performing the non-trivial integration over 
X. Details of this integration can be found in App. The results are 



G'S:^^(^,^o)=21n 



m^rrii 



PI I 



Xq) 



+ 



P.4 



In 



^fxo - Pin - xo) 



2{Pl + m1)[Pl{l-Xo)-mlx,] 



+ 2 



+ 



Pi 



21n(6i)ln 



h^^J \ia{ml + m'i - Pi) 



2 



+ 2^(-l)'=Li2 

k=l 




GIT-* {Pi, Xq) 



2mami{Pl + 2ml 



A, 



-arctan 



+ 



l^ln(xo) - 

ia 

4m^ 



rrii 

mlPi 



Xq] 



m 



^ Ri, 



(Xq) 



ln(66 



2 

a" ia 



{Pi + rr^? 



In 



- Xo 

Xq f (1 - Xq) 



mt 



1 + ^1 
mf 

ml{l-Xo) 



Pi 



Pi + mj 



1 - Xq) 



+ 



P.4 

la 

2Xin 



Xo[Pl-2mlil-xo)] 
2[Pl{l-Xo)-mjxo] 

(4.28) 



with the abbreviations 
bo 



-4mlPl{l-Xo) 



+ Riaixo)]^ + 4ml{Pl + m2)(l - Xo) 



■'2,3 



2mf + Pi ± 



p2 '^0, 

2mfxo-i^i(l-Xo) + 



'M,5 



2{mf -mlT 



bi = 

Ka) 



2m2 - P,2 



/A" 



2m? - Pi + 



Xin 



2ml- Pi T 

XiaRia{Xo){l - Xo) 



Am 



2m2xo - Pl{l- Xq) 



^i?ia(xo)(l - Xq) 



(4.29) 



We have checked that all arguments of the logarithms and dilogarithms in ( [4.28|) lie on the 
first Riemann sheet of the corresponding function for the allowed regions of the various 
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parameters. Although the spin- flip contribution Gia^-\Pia, x) is not IR-singular at x 1, 
we have nevertheless introduced the [...]+ distribution. This is advantageous for a large 
momentum transfer, since in the limit j the complete contribution of Qi^^\Pia, x) 
is contained in the endpoint part G,-™^'*(Pj^, xq). If the limit xq can be set to zero, many 
terms in the endpoint contributions Gf^^\Pl^,Xo) simplify. In order to facilitate the 

application of our results, we list the corresponding results Gj-a"^"* (Pj^) = G|a"^''(Pj^, 0) in 
App. ^ explicitly. 

The final result for the ia contribution to the phase-space integral of the subtraction 
function reads 



a 



ZTT 



X 



dx 



xo 



d<i'0.a(i^l,i^lx) \Mo{pa{x),Mxy,T^f) 
d%,UKL PL 1) Off^liPL X) \Mo{Pa{l),P^{ly, TKf) 



(4.30) 



with //' = ia. As already mentioned above, the phase-space integration over the x- 
dependent momenta has to be performed carefully. First, one has to determine the squared 
CM energy s = {pa{x) + Pb)"^ of the new initial state. The needed scalar product Pa{x)pb 
is obtained upon contracting the first equation of ( |4.19| ) with pi,,^. The product PiaPb, 



which appears on the r.h.s., can be replaced by PiaPb = {fnl + Pia ~ ^ia)/'^^ according to 
the definition of the outgoing momentum Kia- In summary, we obtain s as a function of 
X, PL and Kl, 



.X 



, 2 2., {Pia + '^mlx){ml + Pl 

x{s - mi_ - mi) + " " " 



K. 



2Pi 



{Pl + ml-mi){ml + Pl-K, 



2Pi 



(4.31) 



Thus, it is necessary first to fix x, Pf^, and Kf^ in the phase-space integration over 
d'^Q^ia{K^^, Pf^.x), before the other phase-space-variables can be parametrized, which is 
indicated by the arguments of d$o,ia- 



4.3 Initial-state emitter and final-state spectator 

The case of an initial-state emitter a and a final-state spectator i is kinematically 
identical with the previous one, where the roles played by a and i are interchanged. 
Therefore, we can make use of the variables and auxiliary quantities Xj^, ^m, Pia-, etc. of 
the previous section and adopt the same restrictions on P.^ and Xia- For finite fermion 
masses the auxiliary functions Qai^^^ read 



litl\Pa,Phk) 



iPak)Xi 



2 Xia 



Xia-ml 



Pak 



- 9it'-{Pa,Pi,k), 
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(sub) 



{Pa,Pi, k) 



m„ 



'1 



(4.32) 



In the IR and collinear limits the functions (f4.32|) behave as required in (|2.13| ) and ( p.l4| ). 
The behaviour of Xia and Zia in those hmits is given in (|3.11| ). The auxihary momenta pa 
and pi are constructed as specified in ( [4.17| ), completing the construction prescription for 
the subtraction contribution |A^sub,aiP- 

The separation of the photon phase space also proceeds along the same lines as in 
the previous section, leading to the same kind of convolution over x. In contrast to the 
previous case, in which the singularities appeared for x 1, the collinear singularity 
{Pak —>■ 0) is not restricted to a single point in x, as can be seen in ( p.llj ). Therefore, it 
is necessary to choose the lower limit xq for x small enough so that the complete range of 
X = Xia is covered for small m^; otherwise collinear singularities in (J2x^ — |A^subP) 

remain uncancelled. Since negative values of Xia (if they occur at all) can never lead to 
collinear singularities, our initial restriction xq > is consistent. If effects of 0{ma) are 
consistently neglected, one can simply take Xq = 0, as already done in Section |. 

The analytical integration of the subtraction function over the photonic phase space 
is performed as in the previous section. Hence, we define Gai^r^ according to ( |4.24| ) with 
//' = ai and carry out the simple integration over Zia, yielding 



7 (sub) 



In 



[1 - 2;i(x)][2 - X - Z2{x)] 

/Xia Ria{x) [ 1 - X - ^2(2;)] [2 - X - 2l(x)]_ 

+ Ria{x) {l + x)\n ^ + 



2mix 



2^2 



P2 



1 - 2:2 (x) 



Zi{x) 



7 (sub) 



^(sub)/ p2 



X 



1 — X. 



(4.33) 



For = 0, the function Qa^^\Pia, x) becomes singular at x — » 1. This singularity is split 
off by introducing the [. . .]+ distribution as specified in ( [4.26|) and ([4.27|) with //' = ai, 
thereby defining the endpoint contributions Ga^^\Pia,xo). Actually, this splitting is not 

needed for Gai^-\Pia, x), which is a simple regular function; we proceed this way in order 
to keep the generic description of the polarized and unpolarized cases. The integration 
over X can be performed analytically, yielding 



Gai,+ iPia^ ^0) 



2 In 



m^rrii 



+ 



P.2 



+ 



Pi 



+ In 



2 In 



Xo 



Pt/^-R^aix^ _ 

1^^1(1 -^0) J 4m2(i^i + m?) ^ ^ 
Pj-i^ + 2ml + 7v^ 

1^272 + 2m2xo + l^/\~aRia{Xo) 

_ 272) - 



In 



bn V Ai, 



ln(^ 



+ In 



In 



ia\ 



XoJ 



2ml + 



+ ^(2 + Xo)ln 



Co. 

P2 + 



Pi{l - 272) - 2mlx„ + VA^ft„(l„)J 



mt 



^laRia 



^Pfa ~ V^iaRiaiXo) , 



\ In(co) 



Co 



2E(-i) 

fc=0 



Li2(cfc) 
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^(sub)/p2 N 



with the shorthands 



(4.34) 



Co 



7 



P.2 



ma 



Ci = 61, 



C2,3 



I p2 -r 

2m^ 



C4,5 



-'4,5, 



-P.2 



m/ + ie 



(4.35) 



The variables &o,i,4,5 are defined in ( 4.2S ). Although the variable 7 can become imaginary 
for some values of Pj^, the result for 0^^+ is always real and does not depend on the sign of 

the infinitesimal imaginary part ie. The endpoint contributions G^^^^\Pli) = gI^^^\pI^, 0) 
for the simpler value xq = are explicitly listed in App. |A.2 . 

The result for the ai contribution to the integrated the subtraction function takes the 
same form as in the ia case, but now we have to identify //' = ai in ( [4.30|) . Con- 
cerning the phase-space integration over d^o^ia{Kf^, P^g^,x), the remarks made at the 
end of the previous section apply as well. The IR singularity is contained in the end- 
point part G^^^^{P^^,Xo). However, the coUinear singularity also appears in the function 

0ai^+\Piaj ^o)- Since all singular terms are factorized, the convolution over x itself can be 
carried out with = 0, and with rrif = for light fermions. 

4.4 Initial-state emitter and initial-state spectator 



For an emitter a and a spectator b we keep the definition (|3.21|) of the variables Xab 
and Uab- Moreover, we introduce the abbreviations 



Pab=Pa+Pb-k, S={pa+Pb) 

and the auxiliary function 



Rab{x) 



The function is regular at x — > 1 with RabiX) 
subtraction function we define 



9ab^+{Pa,Pb,k) 



[sx + m^)^ 



4m^ 



9ab!^\pa,Pb,k) 



{Pak)Xab 

ml I 



•^ab 



1 - Xab) 



/^ab 

1 for 

Xab 



Xab = X{s,ml,ml) (4.36) 



(4.37) 



0. For the contribution to the 



2 

Xabma 
Pak 



9^b^-{Pa,Pb,k), 



2{paky 



Xab 



(4.38) 



These functions posses the required asymptotic behaviour in the singular limits, which 
are characterized by ( |3.23| ). The functions Qab^^ are set to zero for Xab < xq, where the 
kinematical lower bound 



Xq> X 



2mamb 



m: 



(4.39) 
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has to be respected. For Xab smaller than x, collinear photon emission cannot occur, and 
the following construction of new momenta would break down. Note that Xq can be set 
to zero in either limit of ^ or — 0, in consistency with our treatment of light 
fermions above. 

As already explained in the case of light fermions in Section |3.3| , the spectator momen- 
tum pf, is kept fixed, but the emitter momentum pa and the momenta kj of all outgoing 
particles other than the photon are changed, resulting in a modification of the total mo- 
mentum Pab = J2j kj. We map the momenta Pa and Pab to 



^ TTTb ['^ ^b 

P^b=P'a+P'b, (4.40) 

leading to the mass-shell relations p^ = and P^^ = P^^. These relations and the 
validity of the required asymptotics in the IR and collinear limits can be checked easily. 
The individual momenta kj are modified by a Lorentz transformation in the same way as 
for light fermions, i.e. we have kj = A^^kj as defined in ( p.26|) . The actual form (|3.27| ) of 
the transformation A'^^ remains valid, but the momenta Pab and Pab of this section have 
to be inserted. This completes the necessary input for the construction of the differential 
subtraction contribution \M.suh,ab\^- 

The separation of the photon phase space is again written in terms of a convolution 
over an auxiliary parameter x, 



J d(j){k,Pab;Pa+Pb)6{x - xo) = J dx J d(j)(^Pab{x);Pa{x) + Pb) J[dk{s,x,yab)]- (4.41) 
The x- dependent momenta 

P^{x)=rAx)+ri (4.42) 

are obtained from pa and Pab upon replacing P^^, by {sx + -|- ml + m^). Thus, they 
coincide with pa and Pab at x = Xab, 

p^^{Xab)=pl P^bi^ab) = P^b- (4.43) 

Note that pl{x) = ml even for x ^ Xab- The measure [dk{s,x,yab)] for the photon phase 
space is derived in App. ^ The result is 

r 1 ry2{x) r2n 

[dk{s,x,yab)] = .3 ^ / ^ ^ dyab / dy.^, (4.44) 

where (f^ is the azimuthal angle of the photon in the CM frame. The integration boundary 
for yab is given by 




yiA^) = ^-^a - xy - (4.45) 
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This relation provides us also with the maximal value Xi of x, for which x = Xab is possible, 



Xi = 1 



(4.46) 



Integrating gah^\pa,Pb, k) M.o [^o,ab', Ti^a] over ip^ results in a trivial factor of 27r. The 



integration of g^j^^!^ over i/ab is also performed easily. Defining 

-2 ry2ix) 



we obtain 

7 (sub) 



Gr^'is,x) 



l + x^ 



xs 



^Vab gai^i\pa,Pb,k), 



(4.47) 



1 — X 



In 



y2[x) 



+ 



2m'ix 



y2[x) yi[x) 



^(sub) 

yab- ' 



S, x) , 



Qab^-\S,X) =1-X. 



(4.48) 



The singularity at x ^ 1, which appears for = 0, is split off by introducing the [. . .] 
distribution. Integration over x yields the endpoint contributions 



s2(l - Xo) 



+ 2(l-xo) + 



In 



m^Kb 



m^s^(l — Xo)^ 



ln((ii 



+ 



-fl — 2xn — X, 



2mi 



2 

(sub) 



— -^o) 



TT 

y 



Gab,- (-5' ^o). 



i(sub) 

where 



1 



G'^b"°^(s,xo) = |(1 -xo)^ 



s + 2ml - \A 



ab 



s + 2ml + 

The integral of the ab part of the subtraction function finally reads 



(4.49) 
(4.50) 



d$l |^sub,ab($l)r = -T^QaCTaQbCTb 



X 



da;^lS^(s,x 



d<l>, 



0,afe(S,Xj 



Mo{pa{x), kj{x); TKa'j 
-J dl>o,ab(s, 1) \Mo{pail),kji^);rKa) 

+ G'i™^^(s,xo) / dl>o,a6(s,l) |A^o(pa(l),^i(l);r«:, 



(4.51) 



where the IR singularity is contained in the function G'^'J,"^''(s, Xq), and the mass singularity, 

which appears in Qai^^{s,xo) as well, is factorized (see Section H). The convolution over 
X itself can be carried out with m^ = 0, and with mf = for light fermions. The relation 
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between the momenta ^Q^ab{x) and the variable x is much simpler than in the case with 



an emitter or a spectator in the final state. Contracting the first equation of (|4.42| ) with 
Pb^^ and cancelling some terms, yields 

s = Pi = -sx + ml + ml + 0{ml;), (4.52) 

i.e. the CM energy of $o,a6(s, is completely determined by the original CM energy 
a/s and X. Knowing s from s and x, it is straightforward to parametrize d^o.afel-s, x). 

5 Applications 

In this section, we compare some numerical results on QED corrections obtained by the 
phase-space slicing method with the ones of the subtraction formalism described in this 
paper. In this context, we mention that we have adjusted all phase-space parametrizations 
to the peaking structure of the integrand in the application of the slicing method. For 
instance, ln(ii^^) is used as integration variable, in order to fiatten the IR pole IZ-E^ 
for small values of the energy of the outgoing photon. Photon emission angles are 
treated in a similar way if collinear photon emission from light fermions can take place. 
These reparametrizations have improved the efficiency of the slicing method considerably, 
whereas such improvements are not necessary for the subtraction method. 

We consider the sample processes 77 — > ff{j), e"7 — > e~7(7), and /x^/^" —>■ PeVein)- 
This choice provides separate applications for the cases zj, ia + ai, and ah of emit- 
ter/spectator pairs //'. 

For the numerical evaluations we take the following set of parameters [ : 

a = 1/137.0359895, 
Mw = 80.41 GeV, Mz = 91.187 GeV, Fz = 2.49 GeV, 

me = 0.51099907 MeV, = 105.658389 MeV, mt = 173.8 GeV. (5.1) 

The weak mixing angle 6^ is fixed by 

cos6'„ = c„ = s^ = Jl-cl- (5.2) 

Mz 

The fermionic couplings to the Z boson are expressed in terms of the vector and axial- 
vector factors 



P f s P 

w,f 

where f = ±1/2 is the third component of the weak isospin of the fermion /. 



5.1 The processes 77 — > ffij) 

The QED and weak corrections to the production of light fermion-anti-fermion pairs 
have been discussed recently in Ref. [ |T^. Details about different variants of phase-space 
slicing and about the dipole formalism presented here can also be found there. In partic- 
ular, the treatment of angular cuts in the phase-space integral is described for the dipole 



formalism. Actually, the subtraction functions of Ref. [ and the ones given in this 
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paper differ by a non-singular factor, leading to a different constant contribution in tlie 



integrated counterparts. We fiave repeated tlie numerics of Ref. [ |T3| for the functions 
defined in tliis paper and found results of the same quality. Using the same number of 
phase-space points in the Monte Carlo integration, which is performed by Vegas [ Q , 
the integration error of the results obtained by phase-space slicing are larger by factors 
of 10-20. 

Here we focus on the case of massive fermions; specifically, we consider the process 

7(fci, Ai) + 7(A;2, A2) ^ tip, a) + t{p, a) [+-f{k, A)], (5.4) 

where ki^2, P, P are the particle momenta, and Ai^25 cr, o' are the corresponding helicities. 
There are two emitter/spectator pairs, tt and tt, both of type ij, and the subtraction 
function ( |2.12| ) is given by 

\Msuh{p,P,k;a,a)\'^ = Qte'^g[^^\p,p,k) \Mo{pi,Pi;t(^,^)\'^ 

+ Qyg^^\p,p,k) \Mo{p2rP2;cr,Ta)\\ (5.5) 

The construction of the auxiliary functions glf^r^ with ij = tt, tt proceeds as described 
in Section [4.1| . In particular, the invariant masses Pij are given by the square of the CM 
energy y/s, 

Ptt = Pit = iP + P + kf = {h + k2f = s. (5.6) 
The pairs of auxiliary momenta {pi,Pi) with / = 1, 2 are obtained from (|4.5| ) upon setting 
Pi = Pi Pj = P and Pi = p, pj = p, respectively. We recall that the spatial parts of the 
spectator momenta and their corresponding auxiliary momenta have the same direction in 
the CM frame, i.e. p||Pi and p||p2- This fact is useful for the implementation of angular 
cuts (see Ref. [0). The integrated counterpart to the differential subtraction function 
( W^ reads 

J d0(p,p, k; ki + k2) \Msnh{p,P, k; a,a)f 

= ^Gi:r\^)[ I dmrPi;h + k2)\Mo{pirPi;ra,a)\' 

+ J d(p{p2,p2;ki + k2)\MoiP2,p2',(^,^^)\'^ , (5-7) 

where we have exploited that the auxiliary functions G'[^'^^''(s) and G[^"^''(s) are already 
fixed by the initial state and coincide. 

For the numerical evaluation of the matrix element Aii of the radiative process 
77 — *• tt7, we apply crossing relations to the result on the related reaction // — > 777, 
which is listed in Ref. [ The phase-space integration is performed by Vegas [ Fi- 
nally, we combine the real-photonic corrections with the virtual photonic corrections, the 



evaluation of which is described in Ref. [ |I5[. The resulting QED correction 5qed to the 



total unpolarized cross section is given in Table |^ for some CM energies y/s. As expected, 
the statistical error of the result of phase-space slicing grows roughly proportional to 
ln{AE/E), where E = yfsjl is the photon beam energy in the CM frame. It is obvious 
that the value !\EjE = 10~^ is still not small enough to guarantee reliable results. For 
smaller values of AE the integration error is again larger by a factor of 10-20 than the 
corresponding error obtained by the application of the dipole subtraction method. 
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v/i/ GeV 


(To/ pb 


Method 


AE/E 


'5qed/% 


360 


0.351422 


Phase-space shcing 


10-2 


1.73343 ±0.00002 


10-4 


1.72936 ±0.00011 


10-6 


1.72932 ±0.00020 


Dipole formahsm 


- 


1.72931 ±0.00001 


500 


0.869434 


Phase-space shcing 


10-2 


0.33592 ±0.00066 


10-4 


0.3309 ±0.0017 


10-6 


0.3307 ±0.0026 


Dipole formahsm 


- 


0.33043 ±0.00013 


1000 


0.428565 


Phase-space shcing 


10-2 


0.1881 ±0.0059 


10-4 


0.184 ±0.013 


10-6 


0.191 ±0.021 


Dipole formalism 


- 


0.17431 ±0.00059 


2000 


0.154450 


Phase-space slicing 


10-2 


0.366 ±0.020 


10-4 


0.346 ±0.045 


10-6 


0.362 ±0.071 


Dipole formalism 




0.3498 ±0.0021 



Table 1: Results on the QED correction (5qed to the unpolarized total cross section of 
77 tt(7). 



5.2 The process e~7 — > e~7(7) 

(i) Moderate scattering energies 

We consider the Compton process 

e"(p,cr) ±7(^7'^) — 



e-(p',a')±7(^'i,A;) [±7(A;^, A'^)], 



where the momenta and helicities are given in parentheses. Owing to the strong polar- 
ization dependence of its polarized cross sections, this process is well-suited to determine 
the degrees of beam polarization of e^ beams. For incoming laser photons and e^ beams 
in the energy region of 1 GeV to 1 TeV, the CM energy is in the MeV range, i.e. the 
CM energy is not large with respect to the electron mass. Details of precision calcula- 
tions, which include the photonic corrections of for such Compton polarimeters 
can be found in Refs. [ |l^, |l^ . In the following we make use of the analytical results on 
the virtual corrections and on the amplitudes for real-photonic bremsstrahlung given in 



Ref. [ IT3] and evaluate the real corrections with the dipole formalism. 

The subtraction function receives contributions of the mixed emitter/spectator types 
ia and ai. Denoting the incoming and outgoing electrons in (|5.8|) by e and e', respectively, 
these contributions are labelled by e'e and ee'. Since both outgoing photons can become 
soft, we have to introduce subtraction functions for each individual final-state photon. 
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Note the IR regions of the two photons are separated in phase space so that the two 
subtraction functions can simply be added. Thus, the full subtraction function reads 

|-MsubP= E \M^kp,p',k[;a,a')\' (5.9) 

1=1,2 

with 

\M!iip,p\ k[- a, = e'githp'.P. \Moipi,p'f, a, Ta')\' 

2 (sub) , 



+ e'g:,Z {P^p',kl) \Mo{pi,p'f,Ta,a')\\ (5.10) 



The functions g^^^^ and g^^^^ are defined in Sections ^]2| and ^73| , where we have to identify 
Pa = Pi Pi = p' 1 and k = k[. The auxiliary electron momenta pi and p[ play the roles of pa 
and pi in ([4.17|) , respectively, where the index / refers to the inserted photon momentum 
k = k[. The subtraction function is completely fixed by the above identifications. Because 
of irta = rrii = m^, we can take Xq = as the lower limit on Xia [see ( [4.14|) ]. 



The integrated counterpart to the subtraction function receives contributions from 



convolutions of the form ( [4.30|) . Owing to Bose symmetry with respect to the interchange 
of the outgoing photons, the two contributions corresponding to the two photons are equal. 
Therefore, we calculate only the integrated subtraction contribution for the photon with 
momentum k'^ and weight this contribution with a factor of 2. Let us first consider the 
phase-space integration in the convolution. The squares of the momenta Pia and Kia are 
given by 

Pi = ip' + K - pf = {k, - k',y = t, Kl = {k, - P^^r = k'i = 0. (5.11) 
Inserting these quantities and rrib = into ( |4.31| ), we obtain 



9 t 2xs + t — 2m: 

5 = ml-- + , , 5.12 

' 2 2/?e'e(x) ^ ' 

for the new squared CM energy used in the convolution over x. The phase-space measure 
d%^ia{Kl, Pi, x) reads 

/d*o.'.(0.'.-) = i(^L/*i |. (5.13) 

The lower limit t^i^{x) on t is determined by two kinematical conditions. Firstly, t cannot 
be lower than — 4£'^, where is the energy of the incoming photon in the original CM 
frame. This condition corresponds to the "edge" of phase space where k[ — >• 0. Secondly, 
the requirement s > ml sets another lower limit on i in the calulation of s from ( |5.12| ) 
for fixed x. Hence, t^i^{x) is the maximum of these two limits. The integration over 
the azimuthal angle ip'2 in the CM frame of pi{x) + k^ yields a factor of 27r owing to the 
rotational invariance of the integrand. In the integrand of the convolution we insert the 
distributions Oi^e^ and ^eeV^ defined in Sections ^]2| and respectively. The endpoint 

contributions G^'e t'* ^^"^ ^^cif^J xo = are taken from App. |A.2| . The auxiliary function 
Piais) is given by 



s — ml 
s — m^ 



(5.14) 
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p p 


(To/ mb 


Method 




Sqeb/% 


+, + 


110.946 


Phase-space shcing 


10-2 


0.4094 ±0.0014 


10-4 


0.4016 ±0.0030 


10-6 


0.4014 ±0.0047 


Dipole formahsm 




0.40131 ±0.00033 


+, - 


65.2608 


Phase-space shcing 


10-2 


0.4996 ±0.0016 


10-4 


0.4921 ±0.0035 


10-6 


0.4898 ±0.0053 


Dipole formahsm 




0.49699 ±0.00092 



Table 2: Results on the 0{a) QED correction 5qeb to the total Born cross section ctq 
of e-7 — > e-7(7) for ^/s — 2.21836 MeV and different degrees of beam polarization Pe 
and P^. 



p p 


(To/ pb 


Method 


AE/E 


A^/rad 




+, + 


90.4372 


IR shcing and 
effective mass factor 


10-2 




5.441 ±0.016 


10-4 




5.416 ±0.031 


10-6 




5.468 ±0.047 


Phase-space shcing 


10-2 


10-2 


5.3783 ±0.0074 


10-4 


10-4 


5.385 ±0.026 




io-<' 


5.151 ±0.055 


Dipole formalism 






5.3588 ±0.0041 


+, - 


12.2425 


IR shcing and 
effective mass factor 


10-2 




15.686 ±0.015 


10-4 




15.685 ±0.026 






15.679 ±0.037 


Phase-space shcing 


10-2 


10-2 


15.655 ±0.0056 


10-4 


10-4 


15.656 ±0.019 


10-6 


10-6 


15.687 ±0.045 


Dipole formalism 






15.649 ±0.011 



Table 3: Results on the 0{a) QED correction Sqed to the total Born cross section (Tq 
of e-7 6-7(7) ^Js = 100 GeV, 20° < < 160°, and different degrees of beam 
polarization P^ and P^. 
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For the squared Born amplitudes A4o(^Pa{x) , Pi{x)^ the invariants s and t correspond 
to the Mandelstam variables s and t as defined in Ref. [ 1^, respectively. 

In Table |^ we give the total cross sections for polarized incoming particles and un- 
polarized outgoing particles for a CM energy that is typical for a Compton polarimeter 
of a future e"'"e~ collider; the beam energies are Eg = 500 GeV and = 2.33 eV. The 
statistical error of the result of phase-space slicing grows with decreasing AE/E^. For 
AE/E^ = 10^^ the infiuence of the finite value of AE is still visible. For the smaller values 
of AE the integration error of the results obtained by the dipole subtraction formalism is 
smaller than the one of the slicing method by at least an order of magnitude. 

(a) High scattering energies 

Compton scattering represents an important reference process in possible future elec- 
tron-photon colliders with CM energies in the GeV to TeV range. In this case the electron 
mass me is small with respect to the CM energy and can be neglected in predictions 
whenever mass singularities are avoided. Detailed discussions of the corresponding lowest- 
order cross sections and the electroweak 0{a) corrections can be found in Refs. [ 0, |T^ . 



In the following we take over the results on the virtual QED corrections given there and 
supplement the calculation of the real-photonic corrections of Ref. [ |n[ by the application 
of the dipole formalism. 

The construction of the subtraction function and its integrated counterpart proceeds 
analogously to the case of finite mg above. One can either expand the above results for 
me — > or make direct use of the general results presented in Section ^ for light fermions. 
There is, however, a difference to the massive case as far as the kinematics is concerned. 
For me — i> 0, exact backward Compton scattering has to be excluded by appropriate cuts 
because of a kinematical u-channel pole in the lowest-order cross section, which is only 
regularized by a finite electron mass. We avoid this singular region by requiring a finite 
angle 6'^ of the outgoing electron with the beam axis in the CM frame. To this end, we 
introduce the step function 

9cnt{0) = 0(0 - ecut)e(180° - ^^eut " &) (5.15) 

and set ^cut = 20° in the numerical evaluation. While the original squared matrix element 
|A^ip is simply multiplied by 5'cut(6'e) in the phase-space integration, the cuts on the 
subtraction function have to be chosen in such a way that the same cuts can be applied 
in the integrated counterpart to the subtraction function. At the same time, one has to 
ensure that the subtraction function still compensates all singularities of \^Al\'^gcut{0e)■ 
Applying the cuts to the polar angles 6'^ ^ of the two momenta p[ in the original CM frame 

fulfills these requirements. Thus, |A^subP (§3) replaced by |A1subp5'cut(^e i)- 

In the limit — > the integrated counterpart to the subtraction function simplifies 
drastically. The boost that relates the CM frames of p + /c^ and pi + = xp + goes 
along the beam axis, and the squared CM energies are related by s = xs. Therefore, the 
phase-space measure and the auxiliary function given above reduce to 

1 ro _ 

d$o,e'e(0,t,x) = / dt, pia{s)=X, (5.16) 

ovrxs J —xs 
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where rotational invariance is already exploited to perform the integration over (^g. Ac- 
cording to the cutting procedure described above, we have to apply the angular cut on 
the angle of p[ in the original CM frame, i.e. we have to transform the polar angle "i^g ^{x) 
of p'l defined in the CM frame of xp + k-y back to the CM frame of p + k^. Denoting the 
angle in the latter frame by 9gi{x), the two angles are related by 

2/ . X X- 1 + (1 +x)cos^?;i(x) 
1 + a; + (x — 1) cos-i?^ 

The cuts are consistently introduced in the convolutions over x if all squared matrix 

elements A^o(pa(0?Pi(0) S^t the factor fi'cut (^0,1(0)5 where ^ is equal to x or 1. 

Let us inspect the IR and mass singularities explicitly. From the results of Ref. [ |l^ 
we deduce that the factor 5qed for the virtual corrections can be decomposed into 
a polarization-independent singular part and a polarization-dependent regular part 

•^QED = -^Ci~t,ml) + 5,U^,X,a',\[), (5.18) 

with the auxiliary function C of ( p.8|) . Since the Mandelstam variable t = {p — p'Y of 
Ref. [ |l^] corresponds to t in the convolution over x described above, the IR and mass 



singularities of the virtual correction exactly cancel against the ones contained in the 
endpoint parts G'^^^I and Gee" + - Therefore, the only uncancelled mass-singular correc- 
tions are the ones contained in the distribution ^ee"+) where they are weighted with the 
splitting function Pff{x) in the convolution over x. 

Table |^ shows our results on the 0{a) QED corrections to the integrated cross section 
for ^/s = 100 GeV and different beam polarizations. The table does not only contain the 
results from the slicing and subtraction methods, but also includes the results obtained 
by a formalism called "IR slicing and effective mass factor". In this approach only the 
IR regions are removed from phase space by cuts, and the collinear poles are regularized 
by applying appropriate factors that replace the poles by the correct mass-dependent 
behaviour. More details about the application of this procedure and of the slicing method 



can be found in Ref. [|T9|. For both slicing variants, the statistical integration errors 
increase with decreasing cut parameters AE/E and A^. Here E is the beam energy in 
the CM frame, and the cut angle A6 defines cones around the electron directions that 
are excluded from phase space. For a cut size of 10~^, the integration errors of the 
different methods are of the same order of magnitude, but at least for the approach with 
effective mass factors the finiteness of the cut is still visible. Therefore, smaller cuts are 
advisable. In this case the superiority of the subtraction formalism becomes obvious. For 
the inspected cuts, there is an improvement of a factor of 2 or more in the integration 
error. PI 



■^We expect that the superiority of the subtraction method is more enhanced if more reahstic cuts are 
appUed. Cutting the electron angle directly, without taking into account a recombination with soft or 
collinear photons in the detector, is a strong idealization. Technically this leads to regions in phase space 
where 5cut(^'c) — 1 ^^nd gcut{d'c ;) = o^' vice versa. For collinear photons these regions shrink to zero, but 
nevertheless induce strong peaks in the integrand. Realistic cuts should avoid such pathologies, leading 
to an improvement in the numerical integration. 
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5.3 The process — > i^ei^eil) 

(i) Moderate scattering energies 

As a final example, we consider the process 



(5.19) 



which is phenomeno logically less important, but — owing to its simplicity — it is well suited 
for demonstrating the application of the dipole formalism in situations with two charged 
fermions in the initial state. At lowest order, there is only an s-channel diagram with 
Z-boson exchange, and the Born amplitude reads 



Mo{p,p',q,q';a, a') 



2{s - M| + iMzFz 



X 



(5.20) 



For our purposes, it is sufficient to describe the Z resonance with the constant experimental 
width Fz given above. The virtual photonic corrections consist of a correction to the 
Znn vertex and the muon wave-function correction. The derivation of the Born cross 
section and the virtual correction is standard and has been performed using the techniques 
described in Ref. [ |^. The results are listed in App. for a finite muon mass m^. The 
bremsstrahlung corrections involve photon emission from the muons in the initial state 
only. The amplitudes for these real corrections have been obtained from the general results 
for /i"^yU~ //7 given in Ref. [|^. 

The subtraction function receives the two contributions g^^-^+ ^ and g^^+]^- ^, which are 
both of type ab, and reads 

\Msuh{p,p',k,q,q';a,a')\'^ = 



2 (sub) 



2 (sub) 



^ip,p',k) \Moip,p,qi,q[;Ta, (r')\' 



9^+J-,rip'^P^f^)\^o{p,p',q2,q2,(^,'r(^')\ ■ (5.21) 



auxiliary variables x = Xab and y = Uab for the g'^li-^J+ 



^. The case of g^^+J- ^ can be treated analogously. The 



In the following, we only describe the construction of the auxiliary momenta p, gi, and 
q[ for the contribution of g'^^^^^ '^'^'^ ^'^^^ ^(''Ub) 

pp' — pk — p'k s 



contribution read 



X 



2ml 



s — 2ml 



y 



pk 
pp' 



(5.22) 



where we have included the relation between x and the two squared CM energies s = 
{p + p'Y and s = {p + p' — kY = {p + p'Y. The subtraction function is consistently set to 
zero for x < Xq with 



Xq 



2ml 

fj, 



(5.23) 



Inserting pa = p, pb = p', and P = Pab = p + p' — k into ( 4.40 ), we get the new momenta 



pp 



4ml{pp'){l 



X 



o^s{x^/Xs + vAi) 



pP = pP+ p'P, 



(5.24) 
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where we have rearranged some terms in order to reveal the behaviour of p in the hmit 



0. The abbreviations and Ag are given by 



A, = A(s, ml, mj) = ^s(s-4m2), A,- = A(S, mj, m^) = -4m2). (5.25) 



For the evaluation of the Born matrix element Aio{p,p', qi, q[) in ( |5.21|) we need a scalar 
product between an initial-state and a final-state momentum, such as p'q[ = PjJ^^q'i ■, in 
addition to the already known quantity s. At this point, the explicit form ( p.27|) of the 
Lorentz transformation A.^^, enters. The calculation of the desired scalar product simply 
requires some contractions among the original momenta and p. 

The evaluation of the integrated counterpart to the subtraction function leads to 
convolutions of the form given in ( [4.51|) . Since the variable x enters the phase space 
^o,ab{s,x) only by the CM energy a/I, the phase-space integration over the squared am- 
plitudes lA^oP is the same as for the lowest-order cross section cro(s) at the CM energy 
Consequently, the convolution can be formulated in terms of lowest-order cross sections. 



and the integrated counterpart to the g^^^-^X r contribution in (|5.21| ) reads 



a 



(sub) 



a (sub) 
+ 2n 



XV As 
s,xo)ao{s,TP-,P+) 



ao(s,rP_,P+) 

(5.26) 



where P:^ are the degrees of beam polarization of the jf^ beams. The factor vA^/vAl 
stems from the flux factors in the transition from squared matrix elements to cross sec- 
tions. 

Table ^ shows some results on the photonic 0{a) corrections to the lowest-order 
cross sections ctq for different /i~-beam polarizations. The /i"*" beam is assumed to be 
unpolarized. The considered CM energy of = 500 MeV is too small for a neglect of 
the muon mass in the non-singular contributions. Therefore, the dependence is 
treated exactly. The results of the dipole subtraction formalism are compared to the ones 
obtained by phase-space slicing, where AE is the cut energy on the outgoing photon, and 
E = A/i/2 denotes the beam energy in the CM frame. Similar to the examples inspected 
previously, for AE/E = 10~^ the influence of the finiteness of AE is still visible at the 
chosen level of accuracy. On the other hand, using the same integration parameters for 
the subtraction method, the improvement in the integration error is between one and two 
orders of magnitude. 

(a) High scattering energies 

Now we turn to high scattering energies and neglect the muon mass whenever possible. 



i.e. we apply the results of Section 3^. In this limit, the virtual correction reduces to the 
simple polarization- independent factor 5qed to the Born cross section do. 



Xvirt 
"QED 



a 

TT 



£(s,m^) 



27r2 



(5.27) 



in agreement with the result given in Ref. [ |2^ on initial-state radiation in e+e — * Z* 



//. Concerning the real correction, the subtraction procedure described above becomes 
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p 


(7o/lO-5pb 


Method 


AE/E 


5qed/% 


+ 


5.74003 


Phase-space shcing 


10-2 


0.3827 ±0.0048 








10-4 


0.352 ±0.011 








10-6 


0.341 ±0.018 






Dipole formahsm 




0.36637 ±0.00022 




9.35981 


Phase-space shcing 


10-2 


0.3275 ±0.0045 








10-4 


0.299 ±0.011 








10-6 


0.292 ±0.017 






Dipole formalism 




0.31238 ±0.00020 



Table 4: Results on the 0{a) QED correction 5qed to the total Born cross section (Tq of 
— > UePeil) for ^/s — 500 Me V and different degrees of //"-beam polarization P_ and 
unpolarized /i'^. 





CTo/pb 


Method 


AE/E 


A^/rad 




+ 


1.32547 


IR slicing and 


10-2 




-4.157 ±0.021 






effective mass factor 


10-4 




-4.331 ±0.055 








10-6 




-4.353 ±0.089 






Phase-space shcing 


10-2 


10-2 


-4.162 ±0.018 








10-4 


10-4 


-4.321 ±0.090 








IQ-*"' 


io-« 


-1.3(3 ±0.22 






Dipole formalism 






-4.29135 ±0.00022 




2.06497 


IR slicing and 


10-2 




-4.168 ±0.021 






effective mass factor 


10-4 




-4.335 ±0.054 








IQ-''^ 




-1.35(3 ±0.087 






Phase-space shcing 


10-2 


10-2 


-4.151 ±0.018 








10-4 


10-4 


-4.257 ±0.091 








10-6 


10-6 


-4.24 ±0.22 






Dipole formalism 






-4.30390 ±0.00020 



Table 5: Results on the 0{a) QED correction 5qed to the total Born cross section (Tq of 
yU+^u- z/cZ/e(7) for ^/s — 50GeV and different degrees of //"-beam polarization P_ and 
unpolarized /i"*". 
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technically simpler, since can be neglected in the kinematics everywhere. Owing to 
the simplicity of the total Born cross section for = 0, 



2 



aoi-s) = (1 - rP+)(l + rP_) ^{v,^ + aj\v^ - ra^y. , ...^ (5.2 



's-Ml + iMzIzI 



the convolution (|5.26| ) over x can be easily performed analytically. Note that the IR and 
mass-singular part of the virtual correction 5qed is again exactly cancelled by the singular 

terms in the endpoint contributions G'['^"^+ , and , . The remaining mass-singular 

contributions are contained in G''^-^J+ + and Q'^^+^J- + and enter the convolution ( |5.26| ) over 
X weighted by the splitting function Pff{x). 

The application of the slicing method additionally requires the analytic treatment of 
photons that are emitted nearly coUinearly from the muon beams, i.e. which have emission 
angles 9 within the ranges 0° < 6* < A6' or 0° < 180° - 9 < AO with A6' -C 1. These 



effects are calculated in the same way as described in Ref. [ 19 1 for initial-state radiation 
in Compton scattering. Details about the variant with effective mass factors can also be 
found there. 

In Table ^ we show the photonic corrections for = 50 GeV, P_ = ±1, and P+ = 0, 
obtained in the small- mass limit 0. The numbers again underline the superiority 

of the subtraction formalism. The integration error is reduced by one to two orders of 
magnitude, without the necessity to look for a plateau in auxiliary parameters, such as 
AE and A9. 



6 Discussion and outlook 



6.1 Features of subtraction methods and the dipole formalism 

As already explained in the introduction, the basic motivation for the development 
of subtraction methods is to avoid singular numerical integrations in the calculation of 
real-photonic (or real-gluonic) corrections. In the previous section, we have compared 
the results for various radiative processes obtained by applying the dipole subtraction 
formalism of this paper with the ones obtained by phase-space slicing. We have found 
that the application of the subtraction formalism typically reduces the integration error 
by an order of magnitude with respect to the results of phase-space slicing, when all 
integrations are performed with the same statistics. As mentioned at the beginning of 
Section ^, the efficiency of the slicing method has been improved by introducing appro- 
priate parametrizations of phase space, whereas such improvements are not needed for 
the subtraction method. 

Moreover, a successful application of the slicing method requires a careful investigation 
of the dependence on the soft-photon cut AE and, if relevant, on the angular cut A9. It 
is necessary to optimize the choice of the cut parameters for all considered observables. 
The integration error roughly grows proportional to the logarithm of a cut parameter if 
the cut becomes small. The optimal choice of cut parameters loosens the cuts as much 
as possible, but still suppresses remnant effects of their finiteness. The optimal set of 
cuts varies with the desired accuracy and, in most cases, also with input parameters, 
such as the scattering energy. In practice, one often tends to choose rather small cuts 
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at the cost of accuracy, in order to be on the safe side. Subtraction methods are not 
plagued by the need of such an optimization procedure. Of course, checking the cutoff 
independence of observables represents a good consistency test of calculations based on 
slicing methods that is not possible for subtraction procedures, but the dipole subtraction 
formalism allows for various other checks, some of which are described in Section |6.3| . 

One of the great advantages of the dipole formalism is certainly its process indepen- 
dence, which distinguishes this approach from most of the other subtraction procedures. 
In this paper, the dipole formalism is worked out for photon radiation in processes in- 
volving charged fermions and any other neutral particles. We stress that all different 
configurations of particle masses and helicities are supported. The subtraction function, 
which removes IR and possible collinear singularities from the differential cross section, 
is constructed in such a way that the transition to the region of small fermion masses 
proceeds smoothly. In other words, there is one subtraction function that interpolates the 
regions of large and small masses. 

Finally, one has to admit that the actual application of subtraction methods, in gen- 
eral, is more involved than the use of phase-space slicing for complicated electroweak pro- 
cesses. The presentation in this paper certainly shows that the application of a subtraction 
procedure can be quite involved for processes with massive particles. The implementation 
of phase-space cuts within subtraction methods is straightforward, but nevertheless can 
be laborious (see also next section). On the other hand, once the procedure is applied 
to a process, such complications are completely overcome, and the advantages described 
above become apparent. 



6.2 Phase-space cuts and distributions 

In the above formulation of the dipole formalism, we mainly concentrated on the 
calculation of total cross sections, but we did not pay particular attention to phase- 
space cuts or to the calculation of distributions. We recall that the difference of the 
differential cross section and the subtraction function is integrated over the full phase 
space $1 of n-|- 1 particles numerically, but the integrated counterpart to the subtraction 
function implicitly contains the integration over the photonic part of phase space, which 
is carried out analytically. The cuts that are applied to the subtraction function have 
to be identical with the ones that are applied and to the integrated counterpart of the 
subtraction function. Otherwise these two contributions will not compensate each other, 
leading to wrong results. In practice, this means that we have to distinguish two types 
of cuts. Firstly, we have the original cuts that are applied to the original differential 
cross section; these cuts concern the full phase space $i of n -|- 1 particles. Secondly, 
we have auxiliary cuts that are applied to the subtraction function and to its integrated 
counterpart; they are defined in the reduced phase spaces $o,//' of n particles. Simple 
examples for the implementation of angular cuts have been described in Section |5.2| for 
Compton scattering at high energies and in Ref. [ r3| for the production of light fermion- 
anti-fermion pairs in photon-photon collisions. 

The calculation of distributions is similar to the application of cuts, since a histogram 
of a distribution is nothing but a series of cuts. Hence, the histogram routine that gen- 
erates the desired distribution during the Monte Carlo integration has to handle each 
column of the histogram in the same way as a cutted contribution to the integrated cross 
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section. Note that this procedure imphes that the original differential cross section and 
the subtraction function may contribute to different columns of the histogram for one 
and the same event. The final result for each column is nevertheless finite, because such 
events are in general far away from the singular regions.0 

6.3 Practical advice 

Subtraction methods offer a number of checks, which are very useful in practice. The 
basic principle of subtraction methods is that all contributions originating from the sub- 
traction function add up to zero in the final result. In the following we describe some 
possible checks for the dipole formalism that are mainly based on this principle. The de- 
scribed checks have been successfully carried out in the applications discussed in Section |^. 

The auxiliary functions gff^l, ^ff'^l^ ^//V '"^^ checked for consistency with- 
out application to a specific process. To this end, one should carry out all integrations 
numerically that have been performed analytically for the derivation of Gffi°l and G^^j}^l. 
Since some of these integrations involve IR singularities, a small photon mass has to 
be consistently used in the numerics. 

For the treatment of specific processes, the construction of the phase spaces $o,//' 
deserves particular care. It can be very useful to compare the corresponding phase-space 
volumes entering the integrations over the phase spaces of n + 1 and n particles. The two 
volumes are obtained as follows: 

(a) in the original integration over d$i we set M.i —>■ 0, Mo — > 1, and gff^l —>■ 1; 

(b) in the integrations of the counterparts over d<l>o,//' we set A4o 1 and use the 
expressions for Off^l and G^^j}'^^ that correspond to gff^l — >■ 1- Those expressions 
can be derived easily, using = 0. 

Note that this phase-space comparison, in particular, represents a non-trivial check on 
the convolutions ([4.30|) in the ia and ai cases, which can be quite complicated for massive 



initial-state fermions. 

In many cases, the phase-space check can be extended by including the full form of the 
functions gfji^l, Sffi^l, and Gj^/'^^, i.e. the only substitutions are A4i —>■ and A4o —>■ 1 
in the phase-space integrations. Owing to the IR and collinear singularities in gff'l, this 
kind of check is not always possible in a simple way. The check is, for instance, useful in 
the ia and ai cases with nia ^ 0. In these cases, all singularities appear for a; — > 1 and 
can be removed by applying the additional cut Xia < 1 — Ax with any small Ax > in 
the integration over d$i. This additional cut has to be incorporated in the convolution 
of Off^l over X, too. The simplest possibility to achieve this is to omit the introduction 
of the [...]+ prescription and to perform the convolution in the range Xq < x < 1 — Ax. 

Of course, many other variants of such consistency checks may be useful in actual 
applications. 



''At the edges of the histogram columns this can also occur for "singular events". The finiteness of 
such contributions is guaranteed by the suppression of phase space for those events. 
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6.4 Generalization to QCD 

In this paper, we have focussed on photon radiation off fermions only. The presented 
formalism can, however, be carried over to gluon radiation for a certain class of processes. 
Consider, for instance, a process that involves a heavy quark-anti-quark pair gg, but 
no other QCD partons. In this case, the gluonic corrections can be obtained from the 
photonic corrections by the replacement Q^a — > 4as/3, and the infinitesimal photon mass 

turns into an infinitesimal gluon mass rrig. Since the IR singularity is abelian, the 
transition to dimensional regularization is performed by the well-known substitution 



where D = A — 2e is the dimension and /i the reference mass of dimensional regularization. 

The results of this paper can also be used to deal with processes involving more than 
one heavy quark-anti-quark pair if the colour flow is treated properly. The colour algebra 
is identical to the corresponding process with massless quarks and can be taken over from 



The presented results do not cover the cases of gluon radiation in which collinear 
singularities are treated within dimensional regularization. This includes real-gluonic cor- 
rections to all processes involving massless partons in the initial state. However, the 
presented results can serve as a starting point for a full generalization of the dipole for- 
malism to QCD with heavy quarks. 

7 Summary 

Following the guideline of Ref. [ ^, where the dipole subtraction formalism is presented 
for NLO QCD corrections involving massless unpolarized partons, we have formulated 
this method for photon radiation off massive fermions. The dipole formalism represents a 
process-independent subtraction procedure that removes all IR and collinear singularities 
from differential cross sections of bremsstrahlung processes. The subtracted singular 
structures are calculated separately, where the integration over the singular regions is 
performed analytically. Consequently, no singular numerical integrations are needed for 
the final result. This advantage distinguishes subtraction formalisms from methods that 
employ phase-space slicing. Slicing methods require a careful optimization of small cuts 
that exclude the singular regions from the numerical phase-space integration. 

Since the consistent inclusion of finite fermion masses turned out to be highly non- 
trivial, we have presented the derivation of the method in a rather detailed way. Our 
formulation, which allows for fermions with definite helicity eigenstates, is applicable to 
processes involving charged fermions and any type of neutral particles. The generalization 
to charged bosons is straightforward. In the limit of small fermion masses, which is of 
particular importance phenomenologically, the dipole formalism simplifies considerably 
and is easy to use. 

In order to illustrate the use of the method in practice, we have applied the dipole 
subtraction method to the processes 77 — >• tt(7), e~7 — >• e~7(7), and — >• i/e^'e(7)- 
Comparing the corresponding results to the ones obtained by slicing methods, we find 




(6.1) 



Ref. [|. 
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improvements in the integration errors of typically an order of magnitude, where Monte 
Carlo integrations are performed with the same statistics in both approaches. 

Finally, we conclude that the dipole subtraction formalism is superior to methods that 
are based on phase-space slicing. Moreover, the presented procedure for photon radiation 
off massive fermions is a first step towards the full generalization of the dipole formalism 
to QCD with heavy quarks. 
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Appendix 

A Special cases 

A.l Light incoming particles 

The case of light fermions in the inital state is of particular importance. For instance, 
it is relevant for e^e~ collisions at high energies, as observed at LEP or the SLC Since 



the ab case with nia^b — ^ is already covered by Section |3.3| , here we concentrate on the 
mixed cases ia and ai with nia — > 0. 

Using the auxiliary parameters and Zia of (|3.10|) , the subtraction functions are 
given by 



9i+{Pa,Pi,k) 

9iL-\Pa,Pi,k) = 0. 
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{Pik)Xia 


.2 


Xia ^ia 








2{p,kyxl 


1 

Zia 


1 




2 


{Pak)Xia 


.2 


Xia 



1 ^in. 



Pik 



9ta!^^il>hPa,k), 



1 - Xi. 



(A.l) 



The lower limit xq on Xia can be set to zero consistently. The construction of the auxiliary 
momenta Pa and pi, which is given in ( [4.17| ) for finite fermion masses, becomes particularly 
simple for nia 0. The result is formally identical with (|3.12| ) for the fully massless case, 
but one should note that p'^ = mf still holds. Using the above relations, the ai and ia 



contributions to the subtraction function |A^subP of (|2.12|) can be constructed easily. 

The calculation of the integrated counterparts to |A^subP also considerably simplifies 
in the limit rua —>■ 0. Since the construction of the auxiliary momenta pi and pa proceeds 
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as in the fully massless case described in Section the convolution over x takes the simple 
form (|3.18|) even for finite mj. The distributions for these convolutions read 
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0. The endpoint contributions are given by 
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(A.4) 



The fully massless limit, in which we additionally have mj — 0, can be read off from the 
above results easily, and we get back the corresponding results of Section ^. 

A.2 Endpoint contributions for a;o = 

In Section^ we have given the endpoint contributions Gjj,'^|(P|j,, Xq) with //' = 
ia, ai, ah for an arbitrary lower limit Xq > and finite fermion masses. In many applica- 
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tions it is possible to set Xq to zero, which simphfies the formulas for these contributions. 
For convenience, we list the results on G^ff%P]f,) = G^^j!'^^{Pff,,0) explicitly. 
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where the abbreviations 6j, Cj, dj, 7 are defined as in Section 

2ml 



2m2 I p2 _ /TT 



(A.5) 

In particular, we have 

(A.6) 



B Derivation of phase-space splittings 



In this appendix, we outline the derivation of the photonic parts of phase space needed 
for the analytical integration of the subtraction function. The emitter/spectator cases zj, 
ia + ai, and ah are kinematically different and are treated separately. 
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Figure 2: Definition of kinematical variables for d(j){pi,pj, k; Pij) in the CM frame of Pij. 

B.l Final-state emitter and final-state spectator 

Tfie photonic phase-space measure [dk{P^j,yij, Zij)] of the ij case is defined in ( [4 .61) . 
In order to derive this measure, we inspect the explicit parametrizations 

J d(p{pi,Pj,k;Pij) = J ^^j j d'^ij J dp- J dp°, 

in the CM frame of Pij. In ( [B.l| ) we have exploited the fact that the momenta pj and pj 
possess the same solid angle Qj in this frame, as a consequence of definition ([4 .51) . The 
phase-space variables of d(f){pi,pj, k; P^j) are illustrated in Fig. ^ The angles 6'^ and ipij 
assign the polar and azimuthal angles of pj with respect to the Pj axis, respectively. The 
particle energies pf and p° can be expressed in terms of the variables Uij and Zjj, 

2mf + P;2(^^ ) 2m] + Pf^il - Vij) 
Pi = ^ 5 Pi = — 1= , (B.2) 

so that ( [4.6|) and (p.lj ) directly lead to [dk{Pfj^yij^ Zij)] as given in (|477| ). The integration 
boundary for the particle energies, which is determined by 



1 > I cos 6ij I 





+ 2p°p° + mj + m| - 


2 


Pil 


|Pi 





(B.3) 



translates into the boundary ( [4.8D of the variables and Zij. 
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B.2 Final-state emitter and initial-state spectator, and vice versa 

The photonic phase-space measure [dk{Pl^,x, Zia)] of the ia and ai cases is defined in 
( [4.18| ). We derive the form of this measure by comparing appropriate parametrizations 
of both sides of ( |4.18| ). On the I.h.s. we factorize the phase space into two two-particle 
phase spaces, 

J d(f){pi, k, Kia, Pa+Pb) = J '^^^2^!- J d(p{pi + k,Kia;Pa+Pb) J d(f){pi,k;pi + k), (B.4) 
and insert the parametrizations 

/ d(P{p, + k,Kia;pa+Pb) = -^^==^== I dPl [dipK, 

I d(j){pi,k;pi + k) = I ^ == I d{p^-paf f d<^^, (B.5) 



where (fx is the azimuthal angle of Kia in the CM frame of pa + Pb, and (f-y is the one of 
the photon in the CM frame of pi + k. Using the relations 

^ -Pi , _ + mf - {pa - pif 

iP^ + ky-pl + mV {p^ + ky-pl + mr ^ 

the integrations over the invariants {pi + k^ and {pi—PaY can be replaced by integrations 
over Xia and Zia. The integration limits -21.2(3^), which are given in ( [4.22[ ), follow from the 
limits on {pa — PiY for fixed {pi + k)^ = {Pia + Pa)"^ = Pfa ~ '"^a ~ Pfal^- They can be 
easily derived in the CM frame oi pi + k. Finally, we use 

in. \ Xj. 



'K{P. + k)\ PL ml) = --—^--^ (B.7) 
on the I.h.s. of ( |4.18| ). On the r.h.s. we make use of the parametrization 

f d(f)(pi{x),K,a]Pa{x) +pb) = f dPfa ^_^= =: j d(pK, (B.8) 

where (fx denotes the azimuthal angle of K^a in the CM frame of Pa{x) + p^. Note that 
s, which is the squared CM energy of Pa{,x) + Ph, depends on s, x, P^ and 

Since we are not dealing with transverse polarizations, but with helicity eigenstates or 

unpolarized configurations, we can assume rotational invariance of M.q ($o,m) around 
the beam axis in the corresponding CM frame oi pa{x) + pb- This implies that the in- 
tegration over the azimuthal angle (px yields a trivial factor of 27r. The integration of 

-A^o (^o,w) over ipx yields a factor of 27r as well because of the one-to-one correspon- 
dence of ifK and (pK^ Inserting the above relations into ( [4. 181) , the photonic phase space 
[dk{PL X, Zia)] can be identified for fixed values of a; = Xia and P^^- The result is given in 

(B). 

^The angles are related by tanipK — f{s, x, Pfa^Kfa) tancp^', if = "^k = is defined in the plane 
spanned by Pa, pb, and Pa- This follows from the fact that components of the direction orthogonal 
to this plane are not affected by the Lorentz transformation that relates the CM frames of pa + Pb 
and Pa + Pb- The integrals over ipx and (px remain unchanged by the transformation: Jg" dipK = 
/p " A(pK\f\/ {cos^ (pK + P sir? (px) = 27r. 
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B.3 Initial-state emitter and initial-state spectator 

The measure [dk{s, x, yab)] for the photon phase space is derived by considering explicit 
representations of the two phase-space volumes d0(. . .) in ( [4.41D . The full phase space of 
Pab and k is parametrized in the CM frame of pa + Pb by 



Ks,Pab,^ 



d<P{k,Pab;Pa+Pb) = ^(1^- J df^7' (B-9) 

where fl^ is the solid angle of the photon. The one-particle phase space of Pab{x) reads 

I d<p{Pabix);paix) +pb) = {2tx)5[pI{x)-PI,) = '^5{x-Xab). (B.IO) 

Putting everything together and expressing the polar angle 9^ of the photon in terms of 
Uab, we get the result (f4.44|) for [dA;(s, x, j/ab)]- The integration boundary ( ^4.45|) on yab is 



determined by | cos^^-y| < 1. 

C Sketch of the calculation of the non-trivial integrals 

In Section ^ we have seen that the analytical integration of the subtraction function 
l-'^subP over the photonic parts [dfc(. . .)] of phase space leads to integrals of a non-trivial 
structure. Therefore, we sketch the calculation of those integrals in this appendix. 

C.l Final-state emitter and final-state spectator 

We first consider the integral for Cf^'!^, as defined in ( |4.9| ), for an emitter i and a 
spectator j in the final state. The integration over the variable Zij is simple and yields 



= r^y^M I 2 In 



l-{l-y)z,{y) 
l~{l-y)z2{y) 

y-[2 + z,{y) + zM][z2{y) - zM] - - zM] 



21/ ' PI, y 



-c^m. (C.l) 



where we have renamed yij to y. The contribution of Gf^^\P'lj) will be calculated below. 
The explicit integral over y involves two types of square roots of quadratic forms in y, 
entering via the limits -21,2(1/) given in ([4.8|) . The limits ?/i,2 are also defined there. Either 
of those roots can be removed by splitting the y range into two pieces: 

(a) ?/i < y < Ay < 1, 

(b) Ay<y<y2. 

The IR singularity is contained in part (a) so that part (b) can be evaluated with rriy = 0, 



replacing the root y — yf hj y. The integration over part (a) is simplified by choosing 
the auxiliary parameter Ay small so that in 0{Ay) the parameter y can be set to zero in 
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the non-singular factors of the integrand. Thus, for 
in this part. Exphcitly we get 
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(C.2) 



The second integrals in both parts are elementary. In the first integrals we remove the 
square roots by the substitutions 



(a) yiVl + X = y + 



yl 



(b) x = y2-y + 
The resulting integrals are of the form 



/ Aw?'- ^ 
(l/2 - y) 2 - y - t/a + 



P.2. 



X2 



dx f{x)\\i[g{x)], 



(C.3) 



where f{x) and g{x) are algebraic functions. Upon decomposing f{x) into partial fractions 
and factorizing g{x), such integrals yield subintegrals that can be expressed in terms of 
logarithms and dilogarithms. A convenient way to obtain compact results is to transform 
the limits xi^2 into and oo in a first step. This is achieved by the substitution ^ = 
{x — Xi)/{x2 — x). The subintegrals that lead to dilogarithms can then be calculated by 
using the standard integral^ 



1 



ln(l + K) 



Z=0,1 



-l)'[Li2(l + + 7]{-au P) ln(l + Pa^)]. 

(C.4) 



Although these steps are straightforward, they nevertheless involve a lot of algebra. There- 
fore, we omit the details. Instead we comment on the IR singularity and the role of the 
parameter Ay. In part (a) the upper limit of the integration over x tends to infinity like 
Ay'^P^j/{mjm'^) for fixed Ay, since the photon mass is infinitesimal. This induces 
terms proportional to \n{Ay /m^) in part (a). On the other hand, part (b) is logarith- 
mically divergent for Ay — » 0. The artifical \n{Ay) terms, of course, cancel in the sum 

^The contributions of the function r]{x, y) — ln(a;y) — ln(a;) — hi{y) are necessary to put the arguments 
of the dilogarithms onto the first Riemann sheet for complex constants ao,i and p. 
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of parts (a) and (b). Finally, we note that we had to exploit some identities for the 
dilogarithms in order to obtain the compact form of the final result ( [4.10|) for G^^^^K 



The calculation of Gj-y"!''' proceeds in a different way. Since this function is IR-finite, 



we can set rriy to zero from the beginning. The defining integral (^4.9|) explicitly reads 



G^'iP^) = JpL /■» d, ^ d. (1^. (C.5) 



Note that only the behaviour of Qij^^ at ?/ ^ is relevant in the IR and coUinear limits [see 
( p.3|) j. Therefore, we have chosen a form of the auxiliary function rij{y) that simplifies 
the integration. We have defined rij{y) in such a way that 



1 - y rjjjy) ^ d_ 

?/2 Rij{y) dy 



Rijiy) 



+ CK). (C.6) 



This choice allows us to perform the integration in ( |C.5| ) over z implicitly upon applying 
integration by parts in the integration over y. The boundary terms of this integration by 
parts vanish, and the result is 

G<->(/^5) = ^ r d, ^ (^^f^f^4(.) - . (C.7) 

'Xii Jo y [ Z2[y) zi{y) J 



where z[2{y) = d-2i,2(l/)/d?/. The final integration over y is elementary. Note that the 
above trick avoids terms such as ln[2;i 2(y)] after the integration over z\ such terms would 
lead to dilogarithms in the final result. 

C.2 Final-state emitter and initial-state spectator, and vice versa 

The integrals for the endpoint contributions defined in ( ^.27| ) for the mixed cases 
//' = ia, ai are calculated in a similar way. Therefore, we outline only the basic steps. 



Inspecting the explicit form of the distributions we find that the integrals ( |4.27| ) 



for Gjj,'^^ again contain two different square roots of quadratic forms in x. Analogously 
to the ij case, we first separate these roots by splitting the range of the x integration as 
follows: 

(a) Xi > X > \ — Ax, with Ax ^ 1, 

(b) 1 — Ax > X > xq. 

Part (a) contains the IR singularity and involves only values of x in the vicinity of 1. 
Thus, in 0{Ax) we can set x to 1 in all non-singular terms of the integral. In particular, 
this replaces the function Ria{x) by 1 + 0{m^) and removes the root implicitly contained 

in Ria{x). In (9(m^) we can replace the explicit root ^ Pf^{l — x)^ — Amfm'^x'^ , which 

appears in zi^2{x) given in ( [4.22| ), by ^ P^^{1 — x)^ — 4m?m^. This root is removed by the 
substitution 



2mjm^ 



ia \ 



Ami ml, 



a;)2 
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The resulting integral is of the form ( C.3 ) and can be reduced to logarithms and diloga- 



rithms, as described above. The IR singularity appears in contributions proportional to 
\n{Ax/m^). In part (b) we can set m-y to zero, since the IR singularity is avoided by the 
finite value of Ax. This eliminates the explicit root in 2:1,2(2^) given in ( [4.22|) . The root 
in Ria{x) is removed by the substitution 

2maA 



(b) y = - -^Ria{x), with A = ^-Pl - mj > 0. 



Note that this substitution is only allowed for A > 0. This is, e.g., fulfilled if P^^ < 0, but 
in general not for all P^^. We evaluate the integrals for the allowed range with A > and 
cover the full parameter space in the final result by analytical continuation in P^^. The 
reduction of the obtained integral, which is again of the general form ( |(J.3| ), to logarithms 
and dilogarithms proceeds as above. However, particular care is needed in the arguments 
of those multivalued functions. As required, the singular ln(Aa;) terms cancel in the sum 
of parts (a) and (b). 

The calculation of Cf^^ is simplified by an appropriate choice of the auxiliary function 
ria{x). Similar to the ij case, we have defined this function in such a way that 



X ria{x) d 



'1 - x)2 RiJx) dx 



X 



+ OK). (C.8) 



Hence, integration by parts can be applied as above, and the resulting integral over x is 
elementary. 

Finally, the calculation of G^^^^^ is trivial. 

C.3 Initial-state emitter and initial-state spectator 

In view of the analytical integrations, the ab case turns out to be the simplest one. 
The integrals of the endpoint parts are given by 

GiC\'^^o)= rdxgi^C\s,x) (C.9) 

Jxo 

with the distributions Gabr^ of ( [4.48[ ). The calculation of G^^^^ is trivial. 

The integral for involves only the square root of the quadratic form in x that 

is contained in the limits 1/1,2(3^) given in ( |4.45| ). Note that the auxiliary function Rab{x) 



does not occur in the integral. As above, we first split the range of the integration over x 
as follows: 

(a) xi > X > \ — Ax, with Ax -C 1, 

(b) 1 — Ax > X > xq. 

In part (a) we can set x to 1 in all non-singular terms, and the root is removed by the 
substitution 



(a) ^ — ^1 + y = l- x + y(l-x)2 

This leads to an integral of the form ( p.3|) , which is evaluated as described above. In part 
(b) we can set to zero, directly resulting in an elementary integral, which is expressed 
in terms of logarithms. 
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D Photonic corrections to — > i^eJ^eil) 



Using the methods described in Ref. [|^, we have calculated the virtual photonic 
corrections of 0{a) for arbitrary muon mass. The one-loop amplitude reads 
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ln{x,) - 1 LMc 



+ 
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w.J(1-75K]}, (D.l) 
where (3 denotes the muon velocity in the CM frame, and Xg is an auxiliary variable, 



P 



13 + 1 



le. 



(D.2) 



The fermion spinors in ( p.l|) carry the same arguments as indicated in the Born amplitude 
M.Q given in (|5.20|) . The spinor chains have been evaluated by applying the Weyl-van der 
Waerden spinor technique, following the formulation of Ref. [ The amplitudes M.i for 
the radiative process /^"'"/i" — > z^e^'e7 are contained in Ref. [ |^ explicitly. The application 
of the slicing method to the real corrections requires the separate calculation of the soft- 
photonic corrections. They are contained in the factor correction d^oit to the Born cross 
section o"o. 
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(D.3) 



which has been deduced from the general results given in Ref. [ ^ . 

The above results can be easily expanded in the limit —>■ 0, which can be used for 
high energies. 
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